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ABSTRACT 



The use of a digital comb filter for the separation of two 
speakers was formulated in previous efforts. An adaptive pitch 
synchronous filter has been developed as an alternative. This 
development stems from the characteristics and structure of the 
speech waveform. In the comb filtering development, there were 
tradeoffs between desired speaker distortion and undesired speaker 
separation. In the development of the adaptive filtering tech- 
niques, the results of these tradeoffs will be examined for 
possible improvements over the comb filtering techniques. 

A series of tests performed on test signal inputs compares 
the performance of both the comb and adaptive filtering systems. 
Included in an appendix are several of the main computer programs 
used in the computer implementation. A series of listening tests 
will be performed with these systems in future work making the 
appendix a necessary part of the thesis. 
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CHAPTER I 
INTRODUCTION 

1.1 Introduction and Goals 

The problem of separating a signal waveform from a noise waveform 
has been the topic of unlimited research for many years. In this 
discussion the problem of interest is one of separating a speech wave- 
form from a "noise waveform" which may be the speech waveform from a 
competing speaker. Two methods for speech enhancement are developed 
in this thesis: comb filtering and adaptive filtering. 

The remainder of Chapter I describes the speech waveform and its 
characteristics in a manner that will be helpful for terminology and 
modeling. In Chapter II a review of the previous systems dealing with 
this specific problem will be covered, and the problem formulation 
introduction for this thesis will be outlined. Chapter III deals with 
the algorithm that was used in determining the fundamental frequency 
of a speech waveform including some of the background research that was 
developed on this method. In Chapter IV the adaptive filtering methods 
are formulated along with a discussion of the rules used and problems 
that should be encountered with the system. Chapter V concentrates on 
the test, signal formulation, processing, and results for the different 
filtering systems implemented. This chapter also deals with the 
comparison methods that will be used between the outputs of the 
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different systems. The computer implementation of the systems will 
be discussed in Chapter VI , and Chapter VII describes the results of 
the actual speech waveforms that were processed. Included in this 
chapter are spectrograms of the inputs and outputs for the systems used. 
The overall results and conclusions are covered in Chapters VIII and IX 
respectively. The Appendices contain computer programs, flowcharts, 
and documentation for the purpose of making the continuing work in this 
problem easier for those persons concerned. 

The goals of this thesis can be summarized in the following 

manner: First, the work that had previously been done was to be 

duplicated. This involved implementing the computer system described 

1 

by Vaden Shields. The purpose of this work was to arrive at a 
starting location before any other work was begun. The system was also 
implemented for the purpose of comparison with any systems that were 
generated by future work. Second, a pitch synchronous adaptive filter- 
ing scheme was developed and implemented on the computer. Now with 
these two completely separate systems, their results could be compared 
to give the third goal. Finally, the ideas of comb filtering and 
adaptive filtering were examined in order to find a limit to their 
effectiveness on the speech enhancement problem. 

1.2 The Speech Waveform and its Characteristics 

Before the speech waveform can be processed by any type of sys- 
tem, the structure, characteristics, and properties of the waveform 
must be understood. A thorough understanding of these items along 
with the definition of several terms is necessary before the various 
processing systems are discussed. 
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Speech production takes place in an area of the body referred to 
as the vocal tract. The vocal tract can be best described as an 
acoustical tube of nonuniform cross-sectional area that originates 
at the vocal cords and terminates at the lips. The vocal tract of 
the average male is 17 centimeters long, and is shorter in the average 
female adult. The nasal tract which lies between the velum and nostrils 
is another nonuniform cross-sectional area tube that may be acoustically 
coupled to the vocal tract by the opening at the velum. 

Within the vocal and nasal tracts lie the major components for 

speech production. The articulators, composed of the lips, jaws, tongue, 

2 

and velum vary with the size and shape of the tract. The velum is 
also referred to as the soft palate and controls the acoustical coupling 
between the vocal and nasal tracts. The vocal cords lie at the lower 
end of the vocal tract. The vocal cords consist of a pair of lips 
made of ligament and muscle, and the opening between the vocal cords is 
referred to as the glottis. The words subglottal and supraglottal 
that are used frequently in the literature refer respectively to the 
areas immmediately below and above the glottal opening. 

In general there are three methods of producing speech waveforms. 
First and most important, voiced sounds are produced by exciting the 
vocal tract with quasi -per iodi c pulses of air that originate in the 
lungs as a steady flow and are chopped into pulses by the vocal cord 
vibrations. Second, the vocal tract is constricted in one area, and 
then by means of forcing the air from the lungs through the constric- 
tion, a turbulent air flow is created. Fricative sounds are produced 

in this manner. Plosive sounds are created by closing the vocal tract 

3 

momentarily, building up a pressure, and then, releasing it. 
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Usual ly sounds that use the vocal cord vibrations for excitation 
are referred to as voiced, and those that do not use the vocal cords 
are referred to as unvoiced. The sounds that are produced with the 
use of the nasal tract are named nasal sounds. For these types of 
sounds the nasal tract is coupled to the vocal tract at the velum. 

The terms pitch, pitch period, and fundamental frequency can all 
be defined at the same time. These terms have been defined in many 
ways by different authors. G. Fant defines them in them in the follow- 
ing manner: 



The basic property of a vocal cord sound is its peri- 
odicity expressed by the duration T Q of a complete 
voice period or by the inverse value of the voice 
fundamental frequency. 

F o * '/T 0 

Fundamental pitch and fundamental frequency are not 
synonymous, but these terms can be used interchange- 
ably due to the close one-to-one correspondence. In 
more strict terminology pitch is a tonal sensation 
and frequency a property of the sound stimulus. The 
duration of a pitch cycle always varies somewhat 
from one period to the next. Such variations are 
systematic determining the intonation in part, acci- 
dental rather than unintentional, but nevertheless of 
importance for the naturalness of human speech. Only 
speaking machines are capable of producing a perfect- 
ly monotonic pitch. 4 



The fundamental frequency of the vocal cords is directly proportional 
to their tension, mass, and the subglottal air pressure. Since most 
adult males have vocal cords that are longer and thicker than adult 
females the fundamental frequency in the male voice is lower than in 
the female voice. 



For a normally speaking voice the fundamental frequency can range 

5 

between 60 hertz and 400 hertz. H. L. Shaffer , through a series of 
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expert men+s, determined that the average fundamental frequency of a 
male speaker is 125 hertz, while the average for a female is 192 hertz. 

Another important term, formant, is defined as the natural fre- 
quency of the vocal tract that corresponds to a resonance or peak of 
energy at various frequencies referred to as resonant frequencies. The 
resonant or formant frequencies depend on the shape of the vocal tract 
and the positions of the articulators. 

There have been several models proposed for the speech production 

6 

apparatus. The model proposed by Oppenheim and Schafer gives a good 
insight into the manner in which the speech waveform is produced. 




MODEL FOR SPEECH PRODUCTION 



7 



Figure 1-1 
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The discrete-time model in figure 1-1 can be examined in several 
sections. The time varying digital filter corresponds to the vocal 
tract, and the coefficients of the filter may be changed to correspond 
to the frequency response of the time varying characteristics of the 
vocal tract. The vocal tract is usually slowly varying in normal speech 
and can be considered fixed for short periods of time (on the order of 
10 milliseconds). The filter can be excited by a train of impulses, 
p(n), or by samples from a random number generator, r(n). In voiced 
speech the impulse train, p(n), is used to excite the filter, and the 
spacing between the samples corresponds to the pitch period. For un- 
voiced speech the sequence r(n) is used as a noise-type source for 

\ 

f i Iter excitation. 

8 

Oppenheim and Schafer describe a factor g(n) that accounts for 
the fact that the actual glottal pulses are not impulses. This addi- 
tional factor can be used with the other two components to describe the 
speech waveform, s v (n), for voiced speech. 

s v (n) = p(n) * v v (n) * g(n) (1.1) 

where: * denotes convolution 

The above expression can also be expressed in its frequency domain 
representation as: 

S v ( e J w ) = P(eJ w ) V v (eJ“) G(eJ“) (1.2) 

As described above the impulse train, p(n), has a spacing between 
pulses Tp corresponding to the pitch period. P(ed t0 ) is also a train of 
impulses separated by intervals of fp = 1/Tp. This shows that the 
energy of the speech waveform will lie in narrow bands centered about 
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the fundamental frequency and its harmonics. 

Whenever the speech is unvoiced, the speech waveform can be des- 
cribed as: 



s u (n) = r(n) * v u (n) * g(n) (1.3) 

Alternatively as before: 

S u (eJ w ) = R(e J w ) V u (ej“) G(eJ“) (1.4) 

There is not much that can be said about the properties of r(n) 
and R(eJ w ) other than it has "noise-like" characteristics. Therefore, 
for unvoiced sounds the spectrum of the signal, s u (n), lacks any type 
of harmonic structure. 

There are certain other aspects of the human speech system that 
should be considered. The ear which is considered to be a part of the 
overall speech system is insensitive to errors in phase of a signal, 
while on the other hand, the ear is extremely sensitive to errors in 
the pitch epochs or pitch periods. Many experiments have shown that 
the pitch of a speech waveform provides its naturalness, and if there 
are any errors in pitch of a processed signal, it becomes immediately 
apparent. 

This discussion has been directed toward the explanation and 
definition of certain aspects of the speech waveform that will be used 
and exploited in the following chapters. It was intended to show how 
the speech waveforms are produced and how these waveforms can be 



mode led. 
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CHAPTER I I 

REVIEW OF THE LITERATURE 



2.1 Shields* System 

Since the foundations of this thesis were taken from another thesis, 
some of the fundamental concepts and ideas used in a thesis by Vaden 
Shields will be covered briefly. Shields in his thesis, "Separation of 
Added Speech Signals By Digital Comb Filtering", suggests that the 
harmonic structure of speech may be the key to the effective removal of 
additive noise. He attacks the problem in the following manner. First, 
he looks into the problem of choosing an optimum unit sample response 
for the filter. Second, he examines the effects of unvoiced speech and 
filter interaction. Finally, he looks into the tracking and detection 
of one signal’s fundamental frequency from a combined signal. The first 
two areas will be covered in this summary. 

The filtering strategy used is simple. Since voiced speech has 

its spectral energy concentrated in narrow harmonic bands, another 

speech signal with the same type of spectrum would have its spectral 

energy concentrated in narrow harmonic bands that would, for the most 

9 

part, not overlap those bands of the first signal. This property 
suggested the use of a comb filter adjusted to allow the first signal 
to pass and the second to be attenuated. 



The comb filter must have several properties when it is used in 
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this application. First, the filter must be time variant. This pro- 
perty evolves from the fact that the pitch period of the speech varies 
with time, and the passbands of the filter must be able to move. 

Second, the filter can only act over a limited portion of the speech 
waveform at one time. This is due to the fact that the approximation 
of periodicity of voiced speech only holds for a short time segment. 

Also, the duration of the impulse response of the filter should be 

10 

finite and short compared to changes in the speech waveform. 

The property above that requires the impulse response of the 
filter be finite implies that the digital filter be nonrecursive or 
FIR. There are several ways to design an FIR filter and the method 
Shields uses is to truncate the infinitely long unit sample response 
with a finite duration data window. If a filter were desired that 
would pass only certain frequencies, the "teeth" or passbands would 
need to be very narrow. One type of filter for this application would 
be a unit sample train in frequency with the spacing of the samples 

fp = 1/MT hertz for - » < f < oo where M is an integer and T is the 

11 

sampling rate. The unit sample response of this filter is a unit 
sample train in time with each sample separated by (M - 1) zeros for 

12 

- °° < n < °°. This filter would be set to pass only harmonics of fp. 

As mentioned above, the infinitely long unit sample response of 
the filter must be truncated in order to create the FIR filter. This 
is done by multiplying the infinitely long unit sample response of the 
filter b.y the unit sample response of the data window. The multipli- 
cation that the above operation implies in the time domain corresponds 
to a convolution of the respective transforms in the frequency domain. 
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The effect of this operation is to create narrow harmonically spaced 
passbands and this is referred to as a digital comb filter. 

The shape of the passbands and the stopbands of the comb filter 

13 

are a function of the type of window that is used. Shields explored 
the effects of four common window functions: the rectangular, the 

Hanning, the Hamming, and the Blackman. 

Another way of looking at what is happening in this process is to 
go back to the step where the filter's infinite duration unit sample 
response and the data window were being multiplied together. Looking 
at these results gives the following equations: 

Suppose x(n) is a finite duration sequence with a Fourier Trans- 
form, X(eJ w ), where: 

X(e Jw ) = Y xCn) e' jwn (2.1) 

n=0 

Now, suppose that x(n) is multiplied by a function f(n), whose unit 
sample response is of the following form: 

m-1 jZini 

f(n) = J_ £ e m ; for - °° < n < 00 (2.2) 

m r=0 

f(n) is a train of unit samples that have (m-1) zeros between each unit 
sample. If z(n) if formed from the product of x(n) and f(n), then, the 
resulting Fourier Transform of z(n) is calculated as follows: 

z(n) = x(n) f(n) (2.3) 

Substituting into equation (2.1) 

N-1 

Z(eJ w ) = l x(n ) f (n) e"J wn 
n=0 



(2.4) 
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Subs+itu+ing equation (2.2) in for f(n) in (2.4) gives: 



N-1 

Z(eJ' w ) = l x(n ) 
n=0 



m-1 j2TTrn 

1 l e m 

m r=0 



r jwn (2.5) 



Performing an interchange of the summations in (2.5) gives: 



m-1 N-1 — j (co - 2ur )n 

Z(e JW ) = 1 l l x(n) e m 

m r=0 n=0 



( 2 . 6 ) 



Equation (2.6) may also be written in the following form: 



Z(eJ u ) = i X{ e j< “ - ^ } 
m r=0 



(2.7) 



Interpreting (2.7) may be easier if considered in the following form: 

Z(e J ’ w ) = 1 fx(e J ' w ) + X(e j(w ~|p) + ... 
m 

+ X(eJ (w • ) \ 

m J (2.8) 

From equation (2.8) it can be seen that if X(ej w ) is periodic on the 

interval {0,2-rr}, then, Z(eJ a) ) is composed of a sequence of teeth that 

are periodic on (0,2 tt_ }, and of the same form as X(eJ w ). 

m 



The difference equation that is implemented by this method can be 
expressed as: 

K 

y(nT) = £ a., x(nT + kmT) 

k=-K 



(2.9) 
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It can be seen that the output is merely a weighted sum of the input 
values separated by mT seconds. The coefficients, a^, are determined 
before the processing starts and are fixed by the input parameter K. 

This dependency only on the value of K may be shown in the following 
manner. 

Using a Hamming Window with length, N = LM, where L = 2K + 1, 
and L = the number of coefficients used, the expression for the window is 

w(n) = 0.54 - 0.46 cos (27rn/N) for 0<_n<N - 1 

= 0 elsewhere. (2. 10) 

Substituting in for the equivalent expression for N: 

w(n ) = 0.54 - 0.46 cos (2ir(iM)/LM) (2.11) 

where n = iM are the only points that have non-zero values. Now, 

w(iM) = 0.54 - 0.46 cos (2ffi/L) (2.12) 

It can be seen that the coefficients have no dependency on the value 
of M or fundamental frequency. This shows how the coefficients may be 
calculated and stored before processing and are always the same after 
the parameter K has been chosen. 

There are several factors or properties of windows that must be 
considered. First, the longer the unit sample response of the window, 
the narrower the passbands of the comb filter will be. This occurs 
when the fundamental frequency decreases. Another consideration is 
the fact that a limit exists for the duration of the window due to 
the changing pitch period. For a fixed window length, the choice of 
window will change the filter characteristics such as; passband width 



14 
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and stopband attenuation. In his thesis, Vaden Shields studied the 
choices for type of window and window length. 

Both of these studies were experimental. The length of the window 
was directly proportional to the parameter K and signals were processed 
with and without an additive noise signal to determine the effects of 
distortion and separation of the original signal. 

The experimental observations agreed with the theoretical results. 

As the value of K was increased (this corresponds to a longer window 
length, more coefficients, and narrower passbands), the separation of 
two speakers improved, but the desired signal became more distorted. 
Shorter window lengths yielded less separation and less distortion. 

Through a series of listening experiments the value K = 3 was chosen 

as a good compromise value from the distortion and separation view- 

15 

points. 

The search for the optimum window type started with the immediate 
elimination of the rectangular window. The rectangular window intro- 
duced significantly more distortion than Hamming, Hanning, or Blackman 
windows. The differences between these three windows were slight, and 

after listening tests the Blackman window was chosen as the best 

16 

overa 1 1 . 

The second area of Shields' thesis dealt with the treatment of un- 
voiced speech. His first approach was to turn off the comb filter upon 
the detection of a pitch period greater than 20 milliseconds. This ap- 
proach seemed to work effectively when used on a signal without any 
noise, but did not effectively separate the desired signal from the noise. 
During unvoiced sections of speech, the noise signal would "pop out" 
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of the processed signal. This effect proved to be distracting to listen- 

17 

ers and caused an impairment to intelligibility. 

Another approach for the treatment of the unvoiced speech seg- 
ments was used. Whenever an unvoiced segment was detected, the comb 
filter would continue to use the last valid pitch period. This ap- 
proach achieved better results for noise suppression, but it tended 
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to distort the desired signal. Shields suggested that other methods 
or approaches for the treatment of unvoiced sections should be investi- 
gated due to the fact that no one method had proven to be superior. 

Shields stated that based on his observations of the results 
this approach to speaker separation works. The performance of the 
system was measured by informal listening tests. The best system 
performance occurred in the recovery of a female speech segment from 
a male - female combined signal. On the other hand, the poorest per- 
formance occurred in the recovery of a male’s speech segment from a 
male - female sum. This can be explained in the following manner. 

This method of comb filtering allows for a variable number of pass- 
bands or "teeth" which depend on the pitch period M. If M is the 
value of the pitch period at a given time, then, the number of passbands 
up to the Nyquist frequency (5 kilohertz) will be M/2. In a female . 
voice the pitch period is generally shorter due to many factors such 
as; size and shape of the vocal tract and art i cu I ators . This shorter 
pitch period would cause fewer passbands to fall within the frequency 
range of the speech spectrum. With fewer passbands present there is 
less chance of falsely recovering an unwanted signal. With the roles 

of the signal reversed, the opposite effect would occur explaining the 
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difficulty in extracting the male's voice from the combination. 
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Durlng the listening tests several factors were noticed that 

caused problems in this system: "unitel I igible unvoiced segments, 

short segments of complete distortion of the signal, and breakthrough 

20 

of the undesired phrase." One cause of the distortion problem was 
attributed to the rapid changes in the pitch period over a short in- 
terval of time. Another cause of the distortion stemmed from the 
inaccuracies in marking the pitch period which was carried out by hand. 

The problem of incomplete or insufficient attenuation of the un- 
wanted speech segment was most bothersome during periods of silence. 
Shields suggested that if periods of silence could be distinguished 

from unvoiced segments, then unvoiced segments and silent segments 
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could be treated separately. 

This approach for comb filtering was also tried with white noise. 

Shields originally thought that if the noise were wideband with most of 

the energy outside the passbands, the si gna I -to-noi se ratio would be 

improved. This method failed to produce the desired result, and it 

changed the wideband noise into noise that was harmonic. This dis- 
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torted the output and produced a signal that had a "reedy" sound. 

Using noise that has highly harmonic, this method showed some 

promising results. The suppression and distortion were much better 

when the fundamental frequency of the noise was below the average 

pitch of the speaker. The "pop ouh" effect could still be observed 

when the passbands of the comb filter moved into the range of the 

23 

frequencies where the noise was present. 

The remainder of Shields’ thesis deals with the automatic extrac- 
tion of pitch period from a combined waveform set by jpstral analysis. 
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Si nee another type of pitch detection will be used in this thesis, this 
section of Shields’ thesis will not be reviewed. 

The conclusions that Shields reached can be summarized in the 
following manner. First, this method will separate the combined 
speech signals of two speakers. Second, this scheme will reduce the 
effects of harmonic noise. Third, possible areas of improvement in- 
clude the treatment of unvoiced segments and periods of silence. The 
listening tests should be conducted in a more sophisticated manner in 
order to reveal errors in the system. Fourth, the area of pitch 
detection could be improved. Also, the choice of coefficients may not 
be optimal, and a better frequency response might be obtained using 
another set of coefficients. Finally, Shields suggests the usage of a 

second comb filter that would use the pitch information of the un- 
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wanted speaker to block or reject the unwanted speaker. 

2.2 Parsons' System 

25 

In a recent paper written by Thomas Parsons, the problem of 
the automatic separation of the simultaneous speech of two talkers is 
approached in a different manner. 

Parsons points out that the brain requires binaural data in 
order to separate the combined signal, and in some cases, a signal is 
received over a channel that does not provide this binaural informa- 
tion. Also, Parsons states that in most schemes, the signal enhance- 
ment exploits characteristics of one signal or statistical differences 

between the signals. The properties or statistics of speech are not 

26 

understood well enough. 
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The method that Parsons uses to separate the two talkers capital- 
izes on the harmonic structure of short segments of speech. In the 
frequency domain the harmonic structure of the speech appears as peaks, 
and the procedure involves removing the peaks of the unwanted speaker, 
and then, taking the inverse transform of the remaining spectrum. This 
is a rather complicated procedure to implement automatically for 
several reasons. The frequency peaks from the two speakers will over- 
lap in some areas. The overlap problem is compounded by the pitch 
variations that occur in natural speech, and this variation causes a 
frequency modulation (FM) distortion of the peaks. Other problems 

occur from the event of the pitch tracks crossing. This may cause 
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errors that lead to jumping between speakers. - 

Parsons uses four rules for detecting an overlap in the peaks, 

and if any of these criteria are met, then, the process of separating 

the peaks begins. By using the a priori knowledge of the peak shapes, 

an estimate of the peak may be subtracted from the combined peaks in 
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order to resolve them. 

The next item addressed is the determination of the peak shapes. 

The largest problem is the effect of the frequency modulation on the 

peak shape. Parsons states that by assuming the pitch variation is 

linear over a short segment, the peaks at the harmonics can be 

approximated by a linear FM ramp. He states that for slow pitch rates, 

the peak shape is approximately the same as in the case of a constant 

pitch period. The frequency modulation of the pitch period is shown 

as a quadrature component with an amplitude proportional to the pitch 

rate, and the peak has a shape that is the second derivative of the 
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in-phase component. This approximation for the peak shape holds for 
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frequencies up +o 3 kilohertz. 

The separation of the two speakers is done by synthesizing the 
spectrum' from a knowledge of the pitch contour and the amplitude and 
phase of every one of the harmonics. The method of synthesis works 
better than the method of subtracting the unwanted speaker's harmonics 
from the combined spectra. Parsons states that if the subtraction 
methods are used, any errors in parameter estimation will result in 
incomplete cancellation of the undesired speaker's harmonics. 

Parsons states that this system is still in its beginning stages, 
and that problems of unvoiced speech still have to be solved. Other 
problems encountered are similar to those mentioned by Shields. These 
problems include pitch detection, tracking one speaker's pitch contour, 
and the areas of speech that have pitch rates much faster than in 
normal speech. 

Both methods summarized in this chapter represent possible solu- 
tions to the speech enhancement problems, and there are probably 
other existing methods. The purpose of this chapter was to acquaint 
the reader with some of the methods that had been proposed, and not 
an attempt to cover all possible methods for speech enhancement. 

2.3 Problem Formulation Introduction 

The problem at this point was to decide how the methods for speech 
enhancement could be improved. The system proposed by Shields was 
used as a guideline because his work had been the only source of 
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information on the subject before Parsons proposed his procedures. 

The primary concern in this thesis was to explore the methods 
used for speech enhancement in voiced sections. The unstructured un- 
voiced speech segments were studied, but because of its lack of struc- 
ture, the procedures were, for the most part, trial and error. Chapter 
IV formulates the method that was proposed as an alternative to the 
comb filter. The method is developed around the structure and time- 
varying properties of the speech waveform. 
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CHAPTER I I I 

PITCH DETECTION PROCEDURES 

The need for a simple and accurate algorithm for determining 
pitch period or fundamental frequency can not be over emphasized. If 
the separation of a speaker and noise is to be done in some small 
amount of time by the methods of comb or adaptive filtering, most of 
the pitch detection algorithms have to be discarded because of their 
complexity and execution time. The accuracy in determining pitch 
period is also very important. The filtering techniques that are 
employed in this thesis use the pitch period information to set the 
time-varying digital filter. These methods of filtering assume that 
the pitch period information is exactly known for a given speech wave- 
form, and therefore, errors in pitch detection downgrade system per- 
formance. An alternative approach for pitch detection is presented in 
this chapter. 

3.1 Research Performed by Henke 

30 

In a recent paper, Henke , discusses the features of an accel- 
erometer signal that measures glottal movements may lead to pitch 
period detection for a speech waveform. The waveforms shown in 
figure 3-1 immediately show why an accelerometer waveform is more 
desirable from the standpoint of pitch epoch marking than the 
corresponding speech waveform. 
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Hen ke states that a small accelerometer can be attached on the 
midline in the suprasternal notch (2 or 3 centimeters below the 
glottis) to measure the varying pressure waveform. He suggested that 
the most prominent and stable feature, the "flyback stroke", be used 
in pitch period detection algorithms. The "flyback stroke", shown In 
figure 3-2, is a feature of the pressure waveform that is character- 
ized by a rapid change from outward to inward acceleration. This 
feature occurs immediately following the maximum outward acceleration. 
The maximum outward acceleration occurs at the instant of closure 
of the glottis or shortly thereafter. Therefore, the negative going 
segment of the waveform that crosses the zero line creates a stable 
point for pitch epoch determination. After these pitch epochs 
have been marked, it is a simple task to determine the pitch period 



or fundamental frequency from these markings. 
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3.2 Methods of Pitch Detection Used 

Henke’s method was used in the work carried out in this thesis. 

A program was written to process a pressure waveform that was ob- 
tained from an acce I erometer at the same time as the speech signal 
was recorded. This was done in a two channel interleaving mode at a 
sampling rate of 10 kilohertz with the speech signal on one channel 
and the acce I erometer signal on the other. 

A program was implemented to find all peaks and zero crossings in 
a segment of the accelerometer waveform. These peaks and zero cross- 
ing locations were stored in two separate arrays, and a threshold was 
used to determine which peaks were the ones that were associated with 
the instance of closure of the vocal cords. The first threshold 
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method consisted of simply setting a constant threshold, by guessing 
before the processing took place. This threshold was used throughout 
the program to determine which peaks were the correct ones. The 
constant threshold approach yielded two types of errors that can be 
described with figure 3-3. First, the problem of false detection or 
detecting a peak that was not associated with glottal closure. Second, 
the false rejection problem had the complementary effect, in that, a 
true peak was not large enough to be detected by the constant threshold. 
The constant threshold method yielded about 40 to 50 percent total 
errors of both types, and the constant threshold was found to be 
speaker and gain dependent. The errors incurred by this method 
suggested a second approach for the threshold determination. 




Figure 3-3 
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The problems in the first approach resulted from the dynamic range 
of the accelerometer waveform. The second approach consisted of set- 
ting a variable threshold. The variable threshold was set by means of 
an energy measurement of the glottal signal over a short time interval. 
The energy, E, was measured in the following manner; 



where N was taken to be 25.6 milliseconds. 

The range of the values for the glottal peaks was divided into 
five levels, and after a careful examination of a few waveforms and 
energy measurements, these ranges were experimentally determined. The 
results of the variable threshold approach were quite encouraging. The 
number of errors made was on the order of 20 percent, and this figure 
could be improved if a more careful statistical study were made to 
determine the optimal threshold settings for a given energy measure- 
ment va I ue. 

The idea of the energy measurement also suggested a concept for 
silence detection. In the thesis done by Shields, he suggests that if 
silence could be detected, then some other approach could be taken to 
improve separation. The silence detector was included in the program 
in order that different schemes of filtering could be tried. A 
silent area was determined when the value of E was equal to zero. 

[The silent areas were marked in the table of pitch period values by 
a minus sign in front of the pitch period valued] 




N- 1 



(3.1) 



n=0 
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Af+er the peaks were detected, zero crossings in the accelero- 
meter signal immediately following the peaks were marked. The zero 
crossing was the termination of the ’’flyback stroke", and as mentioned 
before, is thought to be one of the more prominent and stable features 
in the glottal waveform. These zero crossings were marked by using 
the second channel of the file, since the signal only required one 
channel . 

Since the minimum pitch period is around 2.0 milliseconds, the 
algorithm included a rule that repositioned the pointer in the wave- 
form by 20 samples after a correct peak was detected. This helped to 
speed up the algorithm somewhat. 

There was another problem that is still associated with this 
method of pitch determination. The transitions just before voiced and 
unvoiced sections were not of the same form as during voiced sections 
or unvoiced sections. Even looking at the section of the waveform 
that was considered a transition area did not help in the decision of 
where to put the marks for the pitch period. Each signal had to be 
marked by hand in these areas, and fortunately, there were not many of 
them in a waveform. 

The outputs of the pitch detection included two separate files. 
One file consisted of the pitch marks. This file could be added to a 
speech signal for the filtering of the signal, or it could be added to 
an accelerometer signal to determine whether the pitch detection 
algorithm had made any errors. The second output file consisted of 
pitch period values that were determined from measuring the distances 
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between the marks in the previous file. The file is referred to as 
the pitch table, and it is used in the filtering programs. 

The methods used in this section were intended to make the pitch 
detection as easy as possible. The filtering schemes used in the 
thesis assumed that the pitch information was known exactly. The pri- 
mary function of the automatic pitch marking program was to eliminate 
as many periods as possible that would have to be marked by hand. If 
pitch period detection accuracy had been required in the automatic 
marking program, more improvements would have been made. 
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CHAPTER IV 

ADAPTIVE FILTERING CONCEPTS 



4.1 Development 

The idea for the pitch synchronous adaptive filter was developed 
from the ideas of Siamak Samsam as an alternative to the conventional 
comb filtering techniques. 

The name comb filter probably came about from its frequency 
response characteri sties. The frequency response of a comb filter, 
shown in figure 4-1 consists of a series of "teeth" or "fingers". 




Figure 4-1 
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The implementation used for the comb filter in figure 4-1 was discussed 
in section 2.1. The basic idea behind this implementation was to pad 
the impulse response of a low-pass window function with zeros, using 
the same number of zeros between the coefficients as the pitch period 
of the signal at that particular time. The spacing between the co- 
efficients of the Impulse response was uniform giving the frequency 
response in figure 4-1. 

In this type of comb filtering implementation the length of the 
filter was on the order of five to nine times the length of one pitch 
period, and this is equivalent to a range of filter lengths between 
40 and 140 milliseconds. In a normal speaking voice the pitch period 
or fundamental frequency can change rapidly over several periods 
depending on intonation or stress, and if the filter length is long, 
the comb filter may not handle the correct samples by the procedures 
discussed in section 2.1. The phenomenon discussed above can pro- 
bably be more easily explained by figure 4-2. 




COMB FILTER IMPLEMENTATION 



Figure 4-2 
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Suppose that the signal in figure 4-2 is a speech waveform with 
pitch periods Tpj where p = 0 to 5. It can be seen that the values of 
the pitch periods vary over the length of the filter. Suppose that the 
number of zeros between coefficients has just been changed to cor- 
respond to the pitch period T Q , and all spacings between the samples 
are adjusted so that they correspond to the period T 0 . This results 
in a misplacement of the coefficients in the waveform with respect to 
the beginning of each pitch period. In the case of a its position 
in that particular period is much different from the position of co- 
efficients a Q , a^, and a 2 . 

Now, suppose that the pitch period T , has a smaller error due to 
the pitch detection program. Since all spacings are uniform, the 
error is propagated and compounded in the subsequent spacings. 

With some insight into what the comb filtering was trying to do, 
the adaptive methods were formulated. Suppose that the coefficients 
are rearranged so that the spacings are no longer uniform, but in- 
stead, the spacings will be made equal to the pitch period for that 
particular point in the waveform. Figure 4-3 shows the placement of 
the coefficients for the adaptive filtering methods. For this method 
each coefficient is in the same relative position in its particular 
period. 

This method does not really correspond to a comb filter as shown 
by the frequency response of a sample adaptive filter in figure 4-4. 

The most substantial thing that can be said about the structure of the 
frequency response shown in figure 4-4 is that it has no general 
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Figure 4-3 

structure. The structure of figure 4-4 points out the frequency res- 
ponse of the adaptive filter provides little insight into its 
capabilities in the speech enhancement problem. 




Figure 4-4 
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The output of this system can be expressed in the following 
equation : 



2K 

y(n) = £ aj x( n— I j ) 

i=0 

where: aj = the coefficients or weightings 



I i = 0 

1 2 = m 1 

I j = iri] + m2 



^ I rn^ • m | __ ^ 

and, mj = the pitch period of the i +h period from the 

originating point. 

2K + 1 = the number of coefficients used. 

There is another method of expressing the ideas given earlier. 
This description is aided with the use of figure 4-5. 

The method can be thought of in the following manner. First, 
the waveform is broken into segments according to the pitch epoch or 
positions where the pitch marks occur. Then, these segments are lined 
up as shown in figure 4-5. Then, a weighted average is computed point 
by point as the filter moves in the direction indicated. What this 
operation intuitively does is to develop an average impulse response 
based on the previous several periods. This technique works well in 
conjunction with the assumption that the impulse response of the vocal 
tract is slowly varying. Another important point to be said for this 
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Figure 4-5 

technique is the fact that in the separation problem, the weighted 
samples of the desired speaker will add coherently while the contri- 
butions of the undesired speaker will not add coherently. In the 
comb filtering methods the same statement can not be made as strongly 
when the pitch is changing. Since the averaging process in the comb 
filter looks at different relative positions in each period, the 
reinforcement or coherent addition of the desired speaker is not as 
pronounced as in the adaptive methods. 

It should be pointed out that both methods are identical, and in 
fact, a comb filter when the pitch period is constant. In other 
words, the adaptive filter approaches a comb filter when the 
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pltch period contour is slowly varying and performs a comb filter opera- 
tion when the pitch period contour is constant. When the pitch con- 
tour is fluctuating rapidly, the adaptive filter does not resemble a 
comb filter in the least. 

There is one problem in this method, and this problem with a 
possible solution is discussed in the next section. 

4.2 Overload Problem 

A problem in the adaptive filtering method referred to as the 
"overload problem or phenomenon" can be described with the aid of 
figure 4-6. 




Figure 4-6 
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Suppose that there is a rapid increase in the pitch period in a 
relatively short interval of time. This phenomenon is not rare in 
normal speech. For some speakers a distinguishable feature or 
characteristic is the short pitch period that occurs at the onset of 
voiced speech. This short period can be represented by the segment 
in figure 4-6 labeled as T2. As the filtering moves to the right com- 
puting a point by point average, the method runs into problems at the 
point T2. Up to the point T2 the filtering has occurred in the cor- 
rect manner. After T2 the coefficient 82 begins to use samples in seg- 
ment two for the averaging. This feature disagrees with the concepts 
of the adaptive filtering methods because the coefficients now are not 
in the same relative positions in their respective segments. This 
phenomenon has been named the overload problem and causes an undesirable 
result in the output waveform. The overload phenomemon will be studied 
more in the next chapter. 

There are several proposals for solutions to the problem of over- 
loading. First, the easiest solution from all viewpoints is to ignore 
the problem and let the overloading occur. If the pitch period is 
slowly varying, the prob lem will be mi nor and the coef f i cients will 
not be far away from the same relative position in their respective 
periods. Second, the adaptive filter could be "turned off" at T^ and 
pass the input signal completely for the time interval between T^ and 
T 0 . This idea originated from the fact that the impulse response of 
the vocal tract would be of sufficiently low amplitude in the tails 
of the impulse response so that the filtering operation would not have 
much effect in this area. The final solution proposal, and the one 
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that was used in this thesis was one that padded the short segments in 
the filter's length with zeros in order to make the length of the 
short segments equal to the one in the front of the filter. 

Figure 4-7 explains this procedure. The shortest period, seg- 
ment 3, will have zeros added in the area indicated. 




Figure 4-7 

When the filter reaches the end of segment 4, the same problem occurs 
and zeros are again added to pad out this segment. In effect this is 
the same procedure as "turning off" the individual coefficients asso- 
ciated with the segments that have been padded with zeros. In figure 
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4-7 this would involve setting a 2 equal to zero after T 2 is reached 
and a-j equal to zero after Tj is reached. 

There is one other factor to be considered in this method. The 
gain associated with the filter has been fixed to be a unity gain 
system by: 



2K 

l a k = 1 (4.2) 

k = 0 

When a particular coefficient is set to zero or "turned off", the 
remaining coefficients must be rescaled so that the output is not 
attenuated. 

The third procedure for handling the overload problem makes 
sense from an intuitive standpoint. It says that if a coefficient 
can not be placed in the same relative position in the segment or 
period due to the shortness of that period, then, the contribution of 
that coefficient will only provide erroneous results, and therefore, 
it should not be considered. 

4.3 Rules 

As the filtering algorithms become more complicated and complex, 
the list of rules for handling various situations also grows. Some 
of the more important rules used will be discussed in this section. 

The discussion thus far has dealt only with the procedures for 
handling voiced speech. Normal speech waveforms contain unvoiced 
sections also, and these sections must also be handled. There are 
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three basic areas or situations involved in the processing of unvoiced 
speech sections. These situations are: 

1. The detection of an unvoiced section, and transition 
from voiced speech procedures to unvoiced speech 
procedures. 

2. The procedures that are used when the filter is entirely 
in the unvoiced section. 

3. The detection of a voiced section, and transition from 
unvoiced speech procedures to voiced speech procedures. 

The first situation mentioned above was handled by using the last 
voiced pitch period as shown in figure 4-8. 

In figure 4-8 the speech waveform has been omitted for clarity, 
but the sections are marked as voiced and unvoiced. The values labeled 
rrij are equal to the pitch period of that particular period. In figure 
4-8 (a) the filter position is shown just prior to an unvoiced area 
detection. Figure 4-8 (b) shows the spacing after the first coeffi- 
cient has entered the unvoiced area. In this figure two of the 
spacings have a spacing of m^ while the other two spacings are related 
to their particular periods. In figure 4-8 (c) another spacing 
change is shown so that three spacings are equal to m^ . This pro- 
cedure is repeated until the entire filter lies in the unvoiced area. 

Now the second situation has been entered, and another scheme is 
begun. Since the method of attenuating the input in an unvoiced 
section v/orked satisfactorily in the system proposed by Shields, it 
was used in this system also. In effect, the filtering was terminated. 
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Figure 4-8 



and the input waveform was multiplied by a constant less than one in 
order to produce the output waveform. 



In the third situation the procedure was basically the inverse of 
the method used in the first situation. When a voiced section was 
detected, the filter was initialized to a constant spacing using the 
first period value. As the filter moved further into the voiced 
sections, the spacings were changed to conform to the pitch periods. 



In summarizing this chapter the most prominent feature of the 
adaptive filtering approach is its conformity with the speech waveform. 
The performance of this system will be evaluated in later chapters. 
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CHAPTER V 

TEST SIGNAL SECTION 
5.1 Development of Test Signals 

It was decided that before any new methods for speech enhancement 
were developed, a method of testing the various schemes should first 
be undertaken. Even though psychological listening tests were to be 

conducted in the future, there were just too many variables in the system 
to be handled in a complex listening test. In order to assign some 
types of performance indicies, and in order to fix some of the vari- 
ables, a test input signal was generated and stored for future use. 

The purpose of the test input signal was to alleviate the un- 
certainty of characteristics of the speech waveform, such as, the 
pitch period and the impulse response. When the pitch period and the 
impulse response of a voiced section of speech are known, then, some 
types of input and output waveform comparisons can be made in order to 
determine optimal filter properties. 

A test input signal was formulated in the following manner: 

First, it was decided that enough comparison data could be obtained 
from a relatively short segment of the input signal. This feature 
would provide results that were accurate, and in addition, they could 
be obtained quickly. Second, it was decided that the impulse response 
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of the generation filter would be time-i nvar iant over the duration of 
the signal. This assumption was chosen in order to make the compari- 
son study easier, and it conforms to the properties of speech, in that, 
the impulse response of the vocal cords over short time segments can 
be considered constant. Third, it was decided that the pitch period 
would vary sinusoidally about some mean pitch period. 

The input signal formulation can be described by the following 
diagrams and equations: 

Let x(n) = the test input signal. 

v(n) = the impulse response of the generating filter. 
w(n) = the nonuniformly spaced impulse train corresponding to 
the pitch pulses or the excitation. 

Then, x(n) = v(n) * w(n) or (5.1) 



wfn) w 


v(n) 


x(n) ^ 


^ 




TEST SIGNAL GENERATION SYSTEM 



Figure 5-1 

The impulse response of the generating filter, v(n), was chosen 



to be a damped sinusoid of the following form: 
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where: v(n) = e“ n//ot sin 2irwn 0 < n < L-1 (5.2) 

The model seems to be a very good model for the impulse response of 
the vocal cords. The parameters a and w were variable and were 
chosen so that the input signal that was generated would be a good 
model of speech. 

The length of the impulse response was chosen to be 2T Q , where 
T Q is the value of the mean pitch period. It was felt that in order 
to model speech as closely as possible, the length of the impulse 
response of v(n) should be longer than the pitch periods encountered. 
This would cause some overlap of the impulse response, v(n), in the 
input waveform. 

The nonuni form I y spaced impulse train, w(n), was formulated in 
the following manner: The spacing between the impulses, (the pitch 

period) was chosen to vary sinusoidally about some mean pitch period 
T q . The figure 5-3 shows the procedure for this development. It 
can be seen that this non uniformly spaced impulse train is periodic 
[o, N] . 



on the interval 
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Figure 5-3 

where: 6j = 6 sin 2 tt i for i = 0 to N-l 

IM (5.3) 

Now with these two components the input signal x(n) was gen- 
erated. This gave a waveform whose impulse response, and pitch 
contour were exactly known, and the various methods of speech 
enhancement could be tried on this signal. 

5.2 Processing of Test Signals 

With the test signals that have been described in the previous 
section, several processing methods may now be more closely examined. 
The filtering schemes used in this section were designed so that they 
resembled the actual filtering systems in all respects except for the 
length of the input waveforms used. Since the test input waveform was 
much shorter than the actual speech waveforms, it could be completely 
stored in core memory, and the amount of time needed for the signal 
processing was short. 

There were three filtering systems implemented for the test signal 
input waveform. These systems included the system proposed by Shields, 
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the adaptive filtering systems proposed in Chapter IV with the over- 
load problem, and the adaptive filtering system without the overload 
problem. For the remainder of the thesis the adaptive system with the 
overload problem will be referred to as the adaptive overload system. 

The adaptive system that compensates for this problem will be referred 
to simply as the adaptive system. 

Another problem that received much thought was the problem of com- 
paring these system in some manner to determine the good and bad points 
of both systems. Several tests were decided upon that were related to 
the problems of speech enhancement. In the thesis done by Shields, the 
tradeoffs between desired speaker distortion and undesired speaker 
separation were described. These tradeoffs were also discussed in Chap- 
ter II. It was decided that two separate tests could be conceived for 
this area. 

The first test will be referred to as the input pass test. Figure 
5-4 describes the manner in which this was implemented. By using the 
enhancement system in this manner, the output waveforms should exactly 
resemble the input waveform, or in other words, the system should act 
as an identity system. By examining the inputs and the outputs of the 
systems, the amount of desired speaker distortion could be determined 
for various numbers of coefficients or windows used. 

The second test, described by figure 5-5, will be referred to as 
the input reject test. In this type of filtering the pitch information 
is provided from another input waveform that is totally uncorrelated 
with the test input signal. The idea behind this type of processing 
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Figure 5-4 

is to examine the amount of attenuation introduced into the output 
waveform. Ideally, the output would be zero, because the enhancement 
system would be passing another waveform with a different pitch con- 
tour. This test was designed in order to measure the amount of 
speaker separation that resulted from the various systems. 
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Figure 5-5 
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Ano+her test was devised in order to examine the capability of 
each system when white noise was added to the input signal. Figure 
5-6 shows the implementation of this system. 
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Figure 5-6 

Shields stated that the system proposed in his work seemed to turn 

white noise into a highly harmonic noise in the output waveform that 

32 

was very annoying. This test was designed to determine if the 
adaptive filter had a problem with white noise, and if it did, the 
systems could be compared to determine which one did the best in these 
ci rcumstances. 

Finally, a test could be performed on the sum of two waveforms 
that would duplicate the problems of the overall speech enhancement 
system. It should be pointed out that this test was merely a 
combination of tests one and two, and since the systems were linear, 
superposition would have given the same results. 
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In all cases it was assumed that the entire test signal was 
periodic with a period of twenty pitch periods, and the processing 
only concerned itself with one period of the overall waveform. 

It should be pointed out that the test signals used in these tests 
were not actual speech waveforms. The input waveform was modeled to 
closely resemble the speech waveform during voiced sections. There 
were no comparison tests made on test signals modeled as unvoiced 
speech. In the next section the results of these tests for the differ- 
ent systems implemented will be discussed. 

5.3 Results of Test Signal Processing 

In order that the results of the test signal processing be pre- 
sented as clearly as possible, many figures will be employed in this 
section. There are three basic types of presentation used in this 
section. First, a time domain presentation is used to show how the 
various systems process the input waveforms, and from these figures, 
the input and output waveforms can be compared visually. Second, a 
frequency domain presentation is used. The test signal was formulated 
to simulate a signal that had been sampled at 10 kilohertz. A short 
time spectrum was computed by resampling the input waveform at 5 kilo- 
hertz, multiplying by a Hamming Window, and then, the spectrum was 
computed using a Fast Fourier Transform (FFT) program. The bandwidth 
of the test signal was well below the Nyquist Frequency for the new 
sampling rate. The frequency domain presentation displays logarith- 
mically the square of the magnitude of the Fourier Transform of the 
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signals, and this was computed to eliminate the effects of any phase 
errors that were introduced by the filtering. Finally, an error func- 
tion is shown in some cases, and this error function can be expressed 
in the following equation: 

E = 20 log{|X(e J ’ W ) | 2 } - 20 log{ | Y(e ja) ) | 2 } (5.4) 

where: X(eJ w ) is the Fourier Transform of the input signal 

Y(eJ w ) is the Fourier Transform of the output signal. 

Notice that the error function defined above can be thought of in 
another manner. Consider the linear system shown in figure 5-7 below: 




Clearly, 

y(n) = x(n) * h(n) (5.5) 

Expressing equation (5.5) in the frequency domain: 

Y(e j0J ) = X(e J ’ w ) H(e jw ) (5.6) 

Now, H(e J ’“) = Y(e JU) ) 

X(eJ w ) (5.7) 
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Taking the magnitude and squaring both sides gives: 



H(eJ w ) 


2 = 


Y(eJ' w ) 


2 


= 1 Y(eJ t0 ) [ 2 








X(eJ w ) 




|X(e Jw ) 2 


(5.8) 



Now, taking the logarithm of both sides gives: 

log|H(eJ w ) | 2 = log [ YCe-i 03 ) [ 2 

|X(eJ“)| 2 

log{|Y(ej w )|2} - log{ |X(eJ' w )|2} ( 5 . 9 ) 

Notice that this is almost the same form as the error function. HCeJ 10 ) 
is often referred to as the system transfer function, and the error 
function may be represented as: 



E = log 



1 



|hOT]: 



(5.10) 



Therefore, the error function could be thought of as being related to 
the inverse system transfer function. 

Before proceeding the following discussion will be covered in or- 
der to show how the error function will be used to calculate system 
performance. Consider the system proposed in figure 5-7. Suppose that 
the following result was desired: 

y(n) = a x(n) (5.11) 

Expressing this in the frequency domain and because the system is 
I i near: 



Y(eJ w ) 



a X(eJ w ) 



( 5 . 12 ) 
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Now 



X(eJ'u) 

YCeJ^) 



a 



(5.13) 



loq |x(eJ' w )|2 

\wsm\2 



log 




(5.14) 



E = log{|X(eJ w )| 2 } - log{|Y(eJ w ) | 2 } = logfl] (5.15) 

I a 1 2 



E = log (1) - 2log{|a|} = - 2log{|a|} 



(5.16) 



The above results will be used in later parts of this section for 
comparative analysis between the systems implemented. 

The first test described in section 5.2 was named the input pass 
test. This test was performed with the three basic systems developed: 
Shields', the adaptive overload, and the adaptive. For all systems 
five coefficients were used, and these coefficients were taken from 
the Hamming Window function. 

Figure 5-8 shows the time domain input and output waveforms for 
Shields' system. Several things can be noticed in the output waveform. 
First, there is a definite amplitude modulation throughout the wave- 
form. Second, there are some areas where an "overload bump" occurs 
in the output. Figure 5-9 shows the same presentation for the adap- 
tive overload system. It can be noticed that the amplitude modulation 
does not exist, but there are still some problems with overloading. 

In figure 5-10, the adaptive system results are shown: input and out- 

put waveforms are essentially identical. 



I npu+ 



- 61 - 




Figure 5-8 Input and Output Waveforms - Shield ! s System; Input Pass Test 
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Figure 5-9 Input and Output Waveforms - Adaptive Overload System; Input Pass Test 
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Figure 5-9 (con't.) Input and Output Waveforms - Adaptive Overload System; Input Pass Test 
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Figure 5-10 (con ! +.) Input and Output Waveforms - Adaptive System; Input Pass Test 
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A similar comparison can be made by viewing the spectra of the 
input and output waveforms. Figure 5-11 shows the spectrum of the 
input waveform, while figures 5-12, 5-13, and 5-14 show the output 
spectra of Shields' system, the adaptive overload system, and the 
adaptive system respectively. Figure 5-15 presents the error function 
for each of the systems. 

For this test the system should be performing as an identity sys- 
tem. This says that the value for a_ in equation (5.11) should be 
equal to unity. Substituting this into equation (5.16): 

E = -2 log (1) = 0 (5.17) 

This result is approached in all three systems, but the adaptive system 
comes the closest to matching the ideal error function for this test. 

From all three presentations, the adaptive system clearly performs 
much better than the other two systems on the input pass test. 

The second test performed was named the input reject test. The 
purpose of this test was to determine how well the systems would re- 
ject or attenuate an undesired speaker. This corresponds to the separ- 
ational aspects of speech enhancement. The manner in which this test 
was conducted was simple. A waveform was generated with different 
pitch periods, and the pitch marks from that signal along with its 
pitch table were used to filter the original waveform. 

The second waveform had a pitch period that also varied sinusoid- 
ally about some average value that was identical to the original 
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Figure 5-13 Output Spectrum - Adaptive Overload System; Input Pass Test 
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Figure 5-14 Output Spectrum - Adaptive System; Input Pass Test 
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Figure 5-15 (a) Error Function - Shields' System; Input Pass Test 
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Adaptive System; Input Pass Test 
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waveform. The 6 of the second waveform was changed in order to give 
the different set of pitch periods. It should be pointed out that 
this procedure allows for several of the corresponding periods to have 
the same pitch period. This is synonymous to a crossing of pitch con- 
tours in a pitch period versus time presentation. Ideally, the output 
from the systems should be zero, but this result could not be achieved. 

Figures 5-16, 5-17, and 5-18 show the time domain results of the 
input reject test for Shields', the adaptive overload, and the adaptive 
systems respectively. From these presentations it is extremely diffi- 
cult to establish which system is doing a better job of attenuating 
the input waveform. Each system has certain areas that it does well 
in, and areas that it does not do well in. 

The frequency domain presentation does provide a better picture 
of how the systems performed on this test. The input signal spectrum 
has not been included in this series of figures because it would be 
the same spectrum as the one shown in figure 5-11. The output spectra 
for the Shields', the adaptive overload, and the adaptive system are 
presented in figures 5-19, 5-20, and 5-21 respectively. In all systems 
the general level of attenuation is on the order of five decibels. 

The system shown in figure 5-7 may again be used to measure the 
performance. In equation (5.11), if the value of a_ is less than one, 
this corresponds to an attenuation of the input signal. Consider 
equation (5.16) again. If the value of a_ is less than unity, fhen, 
the logarithm of a_ wi I I be negative making the value of E positive. 

As a approaches zero, the value of E approaches infinity. Therefore, 
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Figure 5-16 Input and Output Waveforms - Shields' System; Input Reject Test 
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Figure 5-16 (con’t. ) I nput and Output Waveforms - Shields' System; Input Reject Test 
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Figure 5-17 Input and Output Waveforms - Adaptive Overload System; Input Reject Test 
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Figure 5-17 (con f t.) Input and Output Waveforms - Adaptive Overload System; Input Reject System 
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Figure 5-18 Input and Output Waveforms - Adaptive System; Input Reject Test 
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Figure 5-18 (con't.) Input and Output Waveforms - Adaptive System; Input Reject Test 
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Figure 5-19 Output Spectrum - Shields' System; Input Reject Test 
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Figure 5-20 Output Spectrum - Adaptive Overload System; Input Reject Test 
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Figure 5-21 Output Spectrum - Adaptive System; Input Reject Test 
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in evaluating the system performance for this test, the magnitude of 
the error function is proportional to the amount of attenuation applied. 
The error functions are given in figure 5-22, and these error functions 
show several characteri sties. The average level of the three error 
functions is approximately the same. Therefore, it may be stated that 
overall the adaptive systems seem to do no worse than the system 
proposed by Shields on the input reject test. 

The next test concerned the performance of the systems in the 
presence of white noise. The white noise was generated by an analog 
noise source, lowpass filtered at 4.7 kilohertz, and then, sampled at 
10 kilohertz. The gain of the noise was selected so that the signal- 
to-noise ratio would be much greater than one. The noise waveform was 
then added to the signal as described in figure 5-6 to produce the in- 
put waveform for the enhancement systems. 

The conclusions from Shields thesis were that the system did not 
perform well in the presence of white noise. The observation that the 
comb filter changed the wideband noise into noise that was highly 
harmonic is exemplified by figure 5-23. The output waveform in figure 
5-23 is quite different from the input signal. The same statement may 
be made concerning the output of the adaptive overload system, which 
is shown in figure 5-24. The output of the adaptive system as shown 
in figure 5-25 is quite different from the outputs of the other two 
systems. The input and output waveforms are almost identical , and 
this is quite encouraging. The adaptive system passes the signal and 
noise combination almost exactly. In the adaptive overload and in 
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Shields' System; Input Reject Test 



87 



N 




ga n i asyvnOs 3am i now 



Figure 5-22 (b) Error Function - Adaptive Overload System; Input Reject Test 
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Adaptive System; Input Reject Test 
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Figure 5-23 Input and Output Wave forms - Sh i e I ds T System; White Noise Test 
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Input and Output Waveforms - Shields* System; White Noise Test 



- 91 - 




Figure 5-24 Input and Output Waveforms - Adaptive Overload System; 
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Figure 5-24 (con't.) Input and Output Waveforms - Adaptive Overload System; White Noise Test 
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Figure 5-25 Input and Output Waveforms - Adaptive System; 
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input and Output Waveforms - Adaptive System; White Noise Test 
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Shields' system, the output waveform is distorted immensely due to 
the white noise. On the other hand, the adaptive system's output is 
noisy, but the distortion of the other two systems is not present. 

Figure 5-26 shows the noiseless input signal and the output of the 
adaptive system with seven coefficients in the filter. It can be seen 
that the effects of increasing the number of coefficients in this case 
i s neg I igi b I e. 

The spectrum of the noise used for this test is shown in figure 
5-27, and the spectrum of the input signal added to the noise is shown 
in figure 5-28. Figure 5-29 is the spectrum of the output from Shields' 
system for the case of the signal with additive white noise. It can 
be seen that the spectrum is composed of distinct harmonic bands that 
were not as prevalent in the other cases. The areas outside the 
harmonic bands have been lost in the filtering process, therefore, 
causing distortion in the output. The output spectrum for the adaptive 
overload system, figure 5-30, resembles the input more closely than 
Shields' system, but there is some energy getting lost as the level of 
the spectrum is less than that of the input's spectrum. The adaptive 
system's output spectrum, figure 5-31, has only minor differences from 
the spectrum of the input. Figure 5-32 shows the error functions of 
the outputs in the same manner that they have been presented in the 
past. 

There is little doubt that the adaptive system performs better 
than the other two systems from the standpoint of this particular test. 
The adaptive system provides an output signal that may be noisy, but 
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Figure 5-26 Input (noiseless) and Output Waveforms - Adaptive System, Seven Coefficients; White Noise Test 
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Input (noiseless) and Output Waveforms - Adaptive System, Seven Coefficients, 
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Figure 5-27 Spectrum of Additive Noise 
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Figure 5-29 Output Spectrum - Shields' System; White Noise Test 
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Figure 5-30 Output Spectrum - Adaptive Overload System; White Noise Test 
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Figure 5-31 Output Spectrum - Adaptive System; White Noise Test 
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Figure 5-32 (a) Error Function - Shields* System; White Noise Test 
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the other two methods provide an output signal that is both noisy and 
d i storted. 

The final test was conducted for the purposes of continuity. 

This test, as pointed out before, is a combination of tests one and 
two. The system proposed by Shields is presented in a time domain 
representation in figure 5-33. The adaptive overload and adaptive 
systems are shown in figures 5-34 and 5-35 respectively. The spectra 
of the three systems are shown in figures 5-36, 5-37, and 5-38 for 
Shields' system, the adaptive overload, and the adaptive systems 
respectively. Upon careful examination of these spectra, the level of 
the spectrum from Shields' system is generally lower at most frequencies 
than either of the adaptive systems or the spectrum of the desired 
speaker shown in figure 5-11. This point seems to agree with the 
phenomenon observed earlier in the input pass test. The amplitude of 
the output from Shields' system in the input pass test was being ampli- 
tude modulated, therefore, the waveform would have to contain less 
energy at some frequencies than an unmodulated waveform of the same 
type. The error functions shown earlier have been omitted for this 
test because they are almost identical and not very benefical for 
compari sons. 

In summary, this series of tests was performed in order to provide 
some other means of system evaluation. It should be stated that the 
fact the adaptive filter did not suppress the noise waveform better 
than the system proposed by Shields is not surprising. The main ad- 
vantage of the adaptive system is fidelity, not rejection. Some of 
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Figure 5-33 Input and Output Waveforms - Shields' System; Combined Test 
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Input and Output Waveforms - Shields' System; Combined Test 
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Figure 5-34 Input and Output Waveforms - Adaptive Overload System; Combined Test 
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Input and Output Waveforms - Adaptive Overload System; Combined Test 
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Figure 5-35 Input and Output Waveforms - Adaptive System; Combined Test 



I nput 



- 112 - 




Input and Output Waveforms - Adaptive System; Combined Test 
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Figure 5-36 Output Spectrum - Shields' System; Combined Test 
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Figure 5-38 Output Spectrum - Adaptive System; Combined Test 



the phenomena and errors observed from these test signals were help- 
ful in implementing the systems that worked with actual speech. When 
the listening tests are performed, a more accurate comparison can pro- 
bably be presented. However, before these tests are performed, other 
than informal listening evaluations and speech spectrograms, these 
results are all that are available for system comparison. 
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CHAPTER VI 

COMPUTER IMPLEMENTATION 
6.1 Description of Computer System 

The computer work for this thesis was done on a PDP-11/45 digital 
computer. This computer is a highly sophisticated and powerful 16-bit 
word machine with 32 K of core memory and dual cartridge disk drives 
capable of holding 1.2 million words each. 

Other peripheral equipment included the Lab Peripheral System, or 
LPS, a storage oscilloscope, a Hewlett Packard 7004 B plotter, and a 
VT05 Alphanumeric Display Terminal. The system was also equiped with 
analog lowpass filters, attenuators, and audio amplifiers. 

The LPS is a modular, real-time subsystem that houses a 12-bit 
analog-to-digital converter (A/D), a programmable rea i time clock, and 
a display contro I I er wh ich includes two 12-bit d i gi ta I -to-ana log con- 
verters (D/A). With the LPS the task of input and output to the com- 
puter was easily handled. 

The LPS real time clock was programmed to sample an input signal 
at 10 kilohertz after a lowpass fi Iter with its cutoff frequency at 
4.7 kilohertz. The computer was set to sample two channels simultan- 
eously, and the speech signal along with the glottal accelerometer 
signal were put onto the disk in the form of a data file. The samples 
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were stored in an interleaving format so that the time reference 
between the two signals would be preserved. The program that conduct- 
ed the sampling was set to sample up to 3.2 seconds of speech with its 
correspond ing accelerometer signal. This amount of time was chosen so 
that most sentences could be accommodated, but with minor modifications, 
the program could be changed to allow any length file. The program 
that conducted the sampling also took advantage of the direct memory 
access (DMA) option of the A/D converter, which allowed the conversions 
to be stored in memory at the maximum rates without processor inter- 
vention. 

The LPS also provided a capable display control. Along with the 
knobs and switches provided, a display program was written that was 
capable of displaying and editing waveforms. This feature was very 
beneficial in correcting pitch markings and in viewing processed and 
unprocessed waveforms. X and Y cursors with knob controls were in- 
cluded with an LED readout of position in order to locate specific 
points in the waveform. 

The routine used for audio output was driven from the D/A convert- 
er in the LPS. The signal from the D/A converter was lowpass filtered 
at 4.7 kilohertz and then amplified. At this point the output signal 
could be recorded on magnetic recording tape for future use. 

6,2 Shields 1 System 

The filtering methods discussed by Shields were implemented on the 
PDP- 11/45 using the standard RT-1 1 Fortran. It was decided that an 
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imp I ernentat ion in fortran would be a slower version, but after algo- 
rithms were perfected, the execution speed could be improved. Flow- 
charts for the systems described in this section are included in the 
appendix. Since the A/D and D/A converters used 12-bits, the input 
and output data values were in integer form although the internal 
multiplications were done in floating point format. 

There were two different systems implemented for the comb filter- 
ing techniques, and the difference in these two methods was the manner 
in which the unvoiced segments were handled. A brief statement about 
these two methods seems necessary before the systems are described. 

The first method, the attenuated input method, stopped the comb fil- 
tering when an unvoiced segment was encountered. The program then 
began to attenuate the input by some specified constant until the 
next voiced section appeared. The second method, the "inertial 
filter", ignored the fact that an unvoiced segment had been encount- 
ered and continued to filter the input using the parameters of the 
last known pitch period value. The name inertial filter resulted 
from the physical aspects of the processing. The filter could be 
thought of as moving through the input samples with some velocity, 
and without any external forces, the filter and its motion would 
remain unchanged as it moved through an unvoiced section. With the 
exception of the manner in which unvoiced speech sections were 
handled, the two systems were basically the same. 

The filter coefficients could be derived from several sources. 

A subroutine was written in order to let the source for the 
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coefficients and the number of coefficients be determined as the pro- 
gram commenced execution. The choices of sources included the four 
window functions that were used in the original system proposed by 
Shields. The window functions were the rectangular, the Hanning, the 
Hamming, and the Blackman. A second choice was the coefficients from 
a lowpass filter designed by the Parks-McClel Ian Algorithm that was 
available on the computer. The final choice was one that allowed the 
user to type in any coefficients desired. 

After the coefficients had been determined, the input buffer was 
initialized with the first section of data. As mentioned in the pre- 
vious section, the waveform had been sampled using two channels. The 
signal was stored on one channel, and the pitch marks were placed on 
the second channel. This type of implementation could be considered 
wasteful from the standpoint of storage, but from the aspects of the 
display, the waveform along with pitch marks could be viev/ed very 
easily. As the filtering commenced the center coefficient checked 
each value on the second channel for the pitch marks, and when a mark 
was encountered, the filter was changed. A pitch table file contained 
the distance between pitch marks and was used to furnish a new pitch 
period value after a mark was encountered. If the value of the pitch 
period was greater than 20 mi lliseconds, the section was treated as 
an unvoiced section. The unvoiced speech procedures described above 
were then implemented depending on the system. 

The pitch table was also marked to denote areas that were silent. 
These si lent areas were detected as mentioned in Chapter I I I from an 
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energy measurement over a segment of data. These silence marks were 
used only to speed up the processing at the beginning and ending of a 
sentence. If a value from the pitch table indicated a silent area, the 
output of the filter was set to zero, and no multiplications were per- 
formed. Procedures for handling the silent areas that appeared inside 
the sentence limits were not thoroughly investigated, and this could 
be a topic for future research. 

From a I I practical aspects and to the best of the author's know- 
ledge, these systems were identical to the ones described by Shields. 
The execution time for these systems was on the order of seventy-five 
times real time. 

6.3 Notch Filters 

Another method that was investigated for the speech enhancement 

problem was a notch filter implementation. Shields suggested that a 

notch filter be used with the pitch marks of the unwanted speaker to 

33 

reject the unwanted speaker. 

In the comb filter implementation, the impulse response of a low- 
pass function was used. Using the principle of duality a notch filter 
implementation yields a frequency response for one of the possible 
filters as shown in figure 6-1. 

The high pass prototype filter was designed by means of the Parks- 
McC lei Ian Algorithm. There were several parameters that had to be 
varied in order to decide what type of prototype fi Iter was needed. 
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Figure 6-1 



These parameters were: 

1 . Fil ter order 

2. Stop Band Width 

3. Transition Width 

The filter order of the prototype high pass filter was limited to 
less than twelfth order. This value was determined from the implemen- 
tation. If a number of'zeros equal to the pitch period were inserted 
between coefficients, then the effective order of the notch filter 
grows quite rapidly. For example, if an eleventh order prototype 
filter were chosen and the pitch period were 100 hertz, then, the 
order of the notch filter would be 1100th order. A filter length of 
this magnitude is on the borderline of the quasi-periodic assumption 
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made for the speech waveform. For the most part, the high pass proto- 
type filter order was between seventh and ninth order. 

The second consideration was the width of the stop band in the 
prototype filter. This parameter can be linked directly to the 
width of the spectrum of the speech waveform at the various harmonics 
of the fundamental. If the stop band is too narrow, then the energy 
at each harmonic from the unwanted speaker will not be sufficiently 
suppressed. On the other hand, if the stop band is too wide, the notch 
filter may suppress the desired speaker and cause distortion. The 
above discussion relates back to the basic tradeoff for speech enhance- 
ment. Again, this tradeoff is the desired speaker distortion versus 
the undesired speaker attenuation. 

The final primary consideration, the transition width, was related 
to the first two considerations already mentioned. The transition 
bandwidth can not be made sufficiently low without increasing the or- 
der of the filter. As the transition width is decreased without chang- 
ing the filter order, the ripples in the pass and stop bands increase. 
Therefore, the limit for the transition width is related to the maxi- 
mum ripple that can be tolerated. 

From the prototype filters that were designed the transition 
widths were relatively large due to the limitations in the filter order. 

After the coefficients from the prototype high pass filter had 
been calculated, these coefficients could be used in the same programs 
as used in the comb filtering implementation. 
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6.4 Adaptive System 

The adaptive filtering systems were implemented much in the same 
manner as those systems mentioned in the previous section. 

Since the methods used to implement the system that corrected for 
the overload problem encompassed the methods used for the system that 
did not, only the former system's description will be covered. 

The initialization was the same as in the previous description, 
but an array was used to hold the information on the spacings be- 
tween the coefficients since this system did not use uniformly spaced 
coefficients. Initially all spacings were set equal to the first 
voiced pitch period. As the filter moved through the input data, and 
pitch marks were detected, the spacing array was continually updated 
to reflect the correct spacing for each coefficient. If an unvoiced 
area were reached, the scheme was to perform the inverse operation 
described above. The coefficient that moved into the unvoiced area 
would retain the last voiced pitch period spacing, and the filter 
continued in this manner until all coefficients were clear of the 
voiced area. In unvoiced areas the input was attenuated as in the 
previous system. 

Another array was used to correct for the overload problem. This 
array contained the initial starting locations in the buffer after the 
filter had been changed. If a coefficient were about to enter an area 
where overloading would occur, that coefficient was set to zero, and 
the rest of the coefficients were rescaled so that the output would 
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have a constant level. The array that contained the starting locations 
was used to detect when an overload situation existed. 

The rest of the system involved only bookkeeping operations that 
were needed to handle the rules used by the system. The system was a 
little slower than the one proposed by Shields and operated on the 
order of ninety times real time. 

In summary it should be pointed out that the fortran implementa- 
tions were not designed for speed. Their basic design philosophy 
was one of user interaction with the programs to ensure correct pro- 
cessing. If speed had been a factor in this implementation, several 
items could have been changed. First, the input waveform could be 
composed of a single channel, and since the data samples were only 
12-bits long, the higher order bits could be used for the pitch mark 
information. With this change in effect longer segments of data could 
be stored in core memory, eliminating several input/output operations. 
The filter implementation could take into account the even symmetry 
of the FIR filter and eliminate one half of the multiplies per output 
point. The addition of assembly language subroutines to do the calcu- 
lational aspects of the program may also speed up the programming. 

With the implementation of these changes, the system begins to lose 
the user interaction capabilities that have proven to be very helpful 



in this work. 
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CHAPTER VI I 

SPEECH WAVEFORM RESULTS 

In order to determine whether or not the systems would perform as 
planned on real speech waveforms, several sentences were processed by 
the various systems. A true evaluation of the performance of the 
three systems involves listening to the processed outputs when speech 
signals are used. These systems will be tested in extensive listening 
tests at a later time, but from informal listening results and spectro- 
graphic analysis will be discussed in this chapter. It should be 
pointed out that it is very difficult to describe the results of the 
various systems with speech waveform inputs. The descriptions and 
comparisons of the systems will be made with characteristic words that 
are not quantitative in the least. 

The initial parameters used in processing the speech waveforms 

originated from the information provided in Shields' thesis. Shields 

concluded that the optimal value for the number of coefficients was 

seven. This number of coefficients provided the best compromise for 

the tradeoff between desired speaker distortion and undesired speaker 

separation. He stated that the Blackman Window provided the best set 

of coefficients from the window functions used, but he also stated 

that the differences in the results between the Blackman and Hamming 

34 

Windows were slight. With these facts in mind, these parameters 
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cou I d remain fixed while the enhancement systems were varied. 

Before the results of the tests similar to those in Chapter V are 
discussed, a few phenomena characteristic to the individual systems 
will be i ntroduced. 

7.1 Comb and Adaptive Filter Results with Speech Waveforms 

There were several phenomena observed in the different systems 
that were not noticed in the test signal section. In the attenuated 
input method of Shields' system, an attenuation constant that was too 
small caused the output waveform to appear chopped. This chopping 
sensation resulted from unvoiced sections that were being attenuated 
by a great amount, and this caused a drastic change in the amplitude 
between a voiced and an unvoiced section. At attenuation constant of 
0.3 seemed to be the smallest value for the constant that allowed the 
chopping to be imperceptible. 

The inertial filter implementation of Shields' system did not pro 
duce the chopping sensation, but a soft reverberation was introduced 
which was audible in quiet. This reverberation could be directly attr 
buted to the fact that unvoiced sections were processed by the comb 
filter with constant parameters. Both the chopping sensation and the 
soft reverberation were phenomena that could be attributed to the 
unvoiced speech segments. Although the time intervals for unvoiced 
sections in normal speech are small, the phenomena occurring in these 
unvoiced sections seemed to carry over into the voiced seel ions and 
disguise the actual results for the voiced sections. 
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The same types of tests performed on the test signals in Chapter V 
were implemented on the actual speech waveforms. These tests were eva- 
luated by informal listening and by spectrograms. Spectrograms will be 
included as a final comparison on the methods that were implemented. 

In these tests the two systems implemented for Shields' method were 
equivalent since the desired speaker consisted of a totally voiced 
waveform. 

For the input pass test all systems performed equally well for the 
most part. The most noticeable characteristic of Shields' system was 
a slight hoarseness in several words. The adaptive systems were char- 
acterized by a slight reverberation. All output waveforms had very 
similar spectrograms with only very minor differences from the spectro- 
gram of the input sentence. 

The input reject test seemed to show that the adaptive systems 
provided about the same attenuation of the undesired speaker, but the 
outputs from the adaptive systems had more reverberation than those 
from Shields' system. 

In the white noise test Shields' system seemed to perform a little 

35 

better. The "reedy" sound mentioned in Shields' thesis was present 
in the output. In the adaptive systems the output had a "buzzing" 
sound that seemed to be more annoying than the distortion in the output 
from the comb filter. Again the spectrograms indicate no major differ- 
ences for this test. 



The combined signal test has been presented with spectrograms, and 
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the combined signal spectrogram is shown in figure 7-1 (a). The two 
waveforms were taken from two different male speakers with approximately 
the same average pitch. The two signals were aligned in time so that 
both speakers would be talking at the same time. This alignment would 
show how well the filtering performed when both speakers were talking. 
Figures 7-1 (b) and (c) show the spectrograms of the two waveforms 
before they were added, and the spectrogram in figure 7-1 (b) has been 
designated as the desired speaker. Figure 7-1 (d) is the spectrogram 
from the output waveform of Shields' system with seven coefficients, 
and a Hamming Window as the lowpass prototype. Finally, figure 7-1 (e) 
is a spectrogram of the adaptive system's output waveform. The output 
waveform in this figure was produced with seven coefficients and a 
Hamming Window as the lowpass prototype also. Both figures 7-1 (d) and 
7-1 (e) are very similar with each one having only minor differences 
from the other. This type of result was normally encountered in other 
sets of sentences that were examined. 

As mentioned earlier in this section, the adaptive systems pro- 
duced a characteristic buzz. The origin of this problem was investi- 
gated by examining the input and output waveforms in the areas where 
the buzzing sound was most prominent. These waveforms are shown in 
figure 7-2. In the output waveform it can be observed that the major 
peak in the input waveform is being clipped and distorted by the 
adaptive system. This phenomenon was not observed in Chapter V when 
working with test signals. In the case of the test signals, the pitch 
period was allowed to change very rapidly in some areas, and the 
adaptive filter performed satisfactorily. It should be pointed out 
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that in the case of the test signals, the impulse response of the 
vocal tract did not change. Therefore, the distortion that is shown 
in figure 7-2 (b) may be due to the increased variation of the impulse 
response in short intervals of time. This phenomenon, fortunately, 
does not occur frequently in normal speech, but it may occur often 
enough to be an annoying problem. 

7.2 Notch Filter Results with Speech Waveforms 

The notch filters were implemented to reject the unwanted speaker. 
Several sentences were processed, and the overall opinion of the notch 
filter was that it did not improve the quality of the output waveforms. 
The notch filter was implemented by itself and in series with a comb 
filter. The notch filter or series combination, at best, performed 
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only as well as the single comb filter by itself. In most cases the 
results of the notch filter or series combination of notch and comb 
were worse than those of the single comb. 

The distortion of the desired speaker was higher in the notch fil- 
ter implementation. This increase in distortion was probably due to 
the wide transistion width that resulted from the design specifications. 
The series combination also had the problem of "harmonic overlap". 

The "harmonic overlap" problem may be explained in the following manner: 
If the teeth of the comb filter are set to pass the desired speaker 
and the notches of the notch filter are set to attenuate the unwanted 
speaker, there may be some areas in the frequency domain when the 
teeth and notches will overlap. This problem can occur even when the 
pitch of the two voices are very different. The overall result of 
this combination of notch and comb has to be lower in performance when 
the "harmonic overlap" occurs. Because of the results from several 
processed sentences the notch filter implementation was not pursued 
further. 

In this chapter the descriptions of the speech enhancement systems 
that were implemented were discussed. This chapter was included for 
the purposes of showing that the systems did in fact provide some type 
of speech enhancement. These systems will be tested in another effort, 
and at that time, some more analytical results may be revealed. A 
demonstration tape for all of the tests described in this section along 
with other sentences was made and may be found in the Digital Signal 
Processing Group Library, M. I. T.. 
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CHAPTER VIII 

RESULTS AND OBSERVATIONS 

The next two chapters will summarize the results and conclusions 
of this thesis. Chapter VIII will deal with the overal I results whi le 
Chapter IX will cover the opinions and general conclusions. 

The general goals of this thesis were fulfilled. The system 
that was described by Shields has been fully implemented on the com- 
puter and was tested thoroughly for correct behavior. As an alterna- 
tive to the comb filtering methods, an adaptive system, has been 
proposed, implemented, and tested. A limited amount of comparison was 
performed with test signals and actual speech waveforms. 

The adaptive methods can be summarized in the following manner: 

The method originated from the time-varying structure of voice pitch. 
The filter was modified to conform to the variations in pitch in or- 
der to pass the desired waveform with as little distortion as possible. 
It was shown in Chapter V that the adaptive system was able to pass 
a time-varying speech- 1 ike waveform when the impulse response of the 
vocal tract was constant. From this viewpoint, it was shown that the 
adaptive methods surpassed the methods formulated by Shields. On the 
other hand, it was stated in Chapter VII that the adaptive systems had 
some problems when the impulse response varied rapidly in short inter- 
vals of time. It should also be pointed out that the adaptive system 
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was more complex and more time-consuming in the computer implementa- 
tion than the system proposed by Shields. 

The adaptive system stressed fidelity of the desired speaker, and 
apparently does no worse on the rejection of the unwanted waveform. 

Finally, it was observed that the notch filter and the notch and 
comb filters operating in series did not provide any substantial im- 
provements in the overall quality of the output waveform. In fact, 
both of these methods showed substantial degradation in some cases. 
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CHAPTER IX 
CONCLUSIONS 

In concluding this discussion the following items should be 
stressed: First, the method proposed in this thesis, the adaptive 

filter, works well in some cases in the speech enhancement problem. 

In other cases, the method has its limitations. The good and bad 
points of the adaptive system were summarized in Chapter VIII. There 
is little doubt that the adaptive filter reduces the level of the un- 
desired speaker in all cases. The desired waveform was sometimes dis- 
torted and occasionally unintelligible. 

Second, the method used for pitch detection as developed in Chap- 
ter III appears to be a very good approach for the pitch detection pro- 
blem. The general pitch detection problem was not considered in this 
thesis, but the entire field of automatic pitch detection algorithms 
has been the topic of much research for many years. This method may 
possibly be improved by first taking the derivative of the glottal 
waveform, and then, using the peak picking algorithm on the resulting 
waveform. There is still one major drawback in this method: The 
transition areas between voiced and unvoiced speech are not easily 
marked whether automatic or manual processing is used. There is a 
general lack of structure in these areas, and a decision has to be 
made as to whether or not a mark is needed and where the mark should 
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be placed. This method of pitch detection also proved to be proficient 
in the automatic detection of si lent areas. Although the maximum bene- 
fit of knowing where these were located was not taken advantage of, the 
silent area detection helped speed up the processing at the beginning 
and ending of the sentences. Overall, this method was generally accur- 
ate and swift. The feasibility of the method may be in question due 
to the fact that the glottal waveform is not as accessible as the 
speech waveform. 

Third, it can be concluded that the notch filter used to reject 
the unwanted speaker does not provide substantial overall improvement 
when both attenuation of undesired signals and faithful reproduction 
of the desired signals are considered. This conclusion is not sur- 
prising, because the notch filter is basically the dual of the comb 
filter. Based on these observations, the notch filter should not 
be used because it introduces the serious distortion of the desired 
waveform. The conclusion is that it seems more feasible to pass the 
desired speaker with as little distortion as possible than to atten- 
uate the unwanted speaker to a higher degree. 

A major observation that should be examined is the question of 
how effective these methods are in speech enhancement. It was con- 
cluded that for normally mixed speech waveforms, there is a limit to 
the amount of enhancement that can be achieved with either of the 
filtering methods discussed. It is believed that even optimally 
designed filters for specific speakers or even for specific sentences, 
would only provide limited improvement because of inherent overlap 
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in the signal spectra. 

The remaining areas open for improvement in this type of approach 
to the problem seem to be in the area of polishing the adaptive system 
so that it approaches this limit. The difficulty in this area lies 
with the fact that much more information concerning the speech wave- 
form will have to be revealed before the techniques can be improved. 
Until that time, alternative methods should be examined for the pro- 
blem of speech enhancement. 
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APPENDIX 



A. 1 I ntroduct i on 

This appendix contains the actual computer programs that were 
used in implementing the systems described in this thesis. They were 
included as an appendix for several reasons. First, if these systems 
need to be implemented again, these programs will provide a reasonable 
starting point. Second, if these systems are to be compared with other 
systems, there will be no question as to whether or not the systems 
described in this thesis have been correctly implemented. Finally, 
these programs were placed in the thesis instead of a separate note- 
book in order that the two sections would not become separated. 

These programs were written in RT-11 Fortran, and all non stand- 
ard fortran routines will be included for completeness. It should also 
be mentioned that some of the features of RT-11 Fortran may not trans- 
late directly into statements usable with another Fortran compiler, 
but for the most part, the RT-11 Fortran is compatible with other 
versions. 

The PDP-11 structure stresses the use of modularity in program- 
ming, and this concept has been used to a great extent in these 
implementations. For each module or program, a brief flowchart is 
included with the computer code. This will probably be a helpful fea- 
ture if the programs in this appendix are to be dissected. 
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Figure A-1 shows a general layout of the overall systems with the 
modules shown as blocks and the computer file names that are used as 
inputs and outputs to these blocks. Figure A-1 is used in the pro- 
cessing of waveforms exactly as shown if there are no pitch errors to 
be corrected. However, if pitch errors occur, then these errors have 
to be corrected as shown in figure A-2. 

These computer programs were not written with speed of execution 
as a primary concern, but with user interaction capability and compre- 
hension for new users as the primary concerns. 
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U K 1 T F < 7 , 7 ) 

( fO^P'aT<* • , * F W I T C t • ( i ) UP - OUTPUT 2 JM EV/£N C-ONpIG.') 

PAUSF ’PAUSE* 

IF ( ISSU'I ( 0 ) . F U . 1 ) CO T.O 3 
CALL ASSIGN (2* ’KU1 jSIGNAL* VYX’ ) 

CALI ASSIGN < 3 »'KK1; SIGNAL* SSS' ) 

CALL ASSIGN (4. •KKIjPITCH.SSS* ) 

GO TO 15 

A URiTF (7.9) 

9 FO*rnT<* * . * NAFiF L| INPUT HUE <F 1 LN A^ • EXT> • ) 

lkjtf:< 7 , i n ) 

1 P FCRPAT ( * » . * * ) 

CALI ASSIGN (2. *KFl jSIGNAL.SSS* ,-14) 

UKITTC 7,H) 

n FOHMflTC *t»p/)MF Ob OUTPUT FILL 1 <FILNaM. LXT>* ) 

URITF ( 7 » 1 U ) 

CALI ASSIGN ( 3 . ,H F 1 ; 5 IGNAL. Oni • , -14 ) 

UKITT (7.12) 

12 f 0 ^ r A T ( * * ,*maFIF OF OUTPUT F T L F 2 < F I L N A N • t X T > * ) 



URITF ( 7,10) 
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CALL ASSIGN < 4 » ’ K Ll ; SIGNAL. 002 • , -14 ) 

]b rLFjfjF F'li.t: ? (2bl,*56,U»IVAR) 

I EFTNL FILL 7 i 2S1 , £ 56 i U . JV aR > 

I LFTNt FILE e ( ?bl ,^b6, U, KVAH ) 
c ape: w thru? 

17 IF ( IVAK.GE.25?) G 0 1C 300 

c ftll the PUFFER. 

ro ?0 )A = 1,10 

F E A P ( ? * I V A R ) ( U A T F Li F (IB + ((IA-3) * ?56 ) ) , IR=1 , ?56 ) 

?0 conjtimuf 

C FILL THE A^^aYS - CHANNEL (4) - UDP POINTS; 

r CHANNEL (14) = EVEN POINTS. 

LC 30 I C = 1 ,2*59,2 
C A T P U F ( IC)=DATdllL ( 1 C ) 

CAlnUF ( IC + 1 ) =0 

30 CONTINUE 

FO 40 I D= 1 , 2 559 , 2 

if ( issui( n .Eu.o) [ r tbup ( j n ) = datbhf<id + i) 

IF ( ISSWIM ) ,cu.0> L A T 13 U F ( ID + 1) = 0 
lF(ISS^I(i j.rQ. 1 ) I / TPUF ( I U ) = o 

IP (ISS w I ( i)..EU.l) LATBUF(1D + 1) = DA 1 BUF ( I D+ 1 ) 

4 (' CONTINUE 

c L'RTTF. THr OUlPUT ONTO T FKl SPECIFIED FILES, 

I 0 P 0 IE = 1,10 

• i,filTF.' ( 3 • J\/AR ) ( GATBlE ( IF + ((IF-1) * 256 ) ) , JF = 1 , 25*> ) 



f ONTIf'UF 
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io r-n xr= 1,10 

UKlTr«H'KV/\R) (RATLIFdH + ((IP-1) * ?b6 ) ) « IH = 1 , 2?b ) 
MJ cOMTTNUf 

GO TO 17 

600 rKOriLE 2 ! Tl'RIM l*']f CS OFF. 

fNOFILE 3 
t l\l 3 F I L. L 4 
Si OP 

f r j d 
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PIT V02 




d> 



FIGURE A- 5 



PIT V02 (CONT.) 
© © 




FIGURE A- 5 (CONT.) 
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c; 

( 

c 

C llUr: PIT v/°2 7b[220 

r 

C SURROIT UFS USED: 

C PEKS 

c zeros 

c 



C PARAMETERS USED a m D PFF lllTiOM; 

c jFlag: this is u<?eu in denote am overflow im tht buffek. 

IT UCCIIKS WMFN A PEAK IS MULCTED AT THE END OF A 
RUFFf.R, aPD THE ( CRKESPONPIUO ZERO CROSSING 
HAS MOT BFEN F Olif [' , WHEN jF e AG IS SET (=1) 

THE PROGRAM I OOKS FOR A ZERO CROSSING IMMEDIATELY 
AFTER a RiiFFTK S| TFT HAS BEt M PERFORMED. 



c zflac: this flag is usE[ to denote a silent area 



c: 

c 

c 



and l HEN SET (=1) A PITCH PERIOD IM THE PITCH 
TABLE JS GIVEN t NEGATIVE SIGN TO DENOTE SILENCE, 



C T E ] riALl7ATlON 
C 



I NT rC-ER DA I BUF ( ?EG(| ) , A( 12B0 ) . B( 3 ??0 ) , 3UF( 2S60 ) 
• INTEGER SUM( S) i ZFLAC . STaR'» TH, E N D 1 



COMMON DATtUlF 
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r at a b jf , i , j, s i ahi go *o « i . 1 » 3 / 

C A I A I V/'.R , JVaR , KV f\(\ , IR, ZFL AG, Jf L A G/ 2 , 2 » 2 t 1 « 0 « 0 / 

C A LI ASSIGN (2« *KK1 ;PTTCH.SSSt ) 

CALL ASSIGN ( 3 , • KK 1 . | J TCH. I AR • ) 

call assign ( < 4 , • kki jpitch. uat « > 

fEFTNt FlLt 2 (2bl ,15G,U, IUAP) 

LEFlNE FlLt" ? (1000, 2, U, JVAR) 
rtFTME FILE u ( ?bl # F 56 » U , K V AR ) 

30 IF ( T VAp , Ew. ?52 ) Gc, TO 300 ! Po WF MUlT? 

1=3 IMP, RtSLT THE POINTER. 

L r ILL THE PUFFER 

to ? n I I = 1 , 1 0 

R L A P ( 2 * jVAp) ( U A 1 PUj ( J + ((IT-1) * 25& ) ) * J = 1 ♦ 25ft > 

2 P CONTINUE 

C FIND THE PEAKS aNP ZERO CROSSINGS. 

CAL) PEKS ( a,Ry IVAH,SUi v !) 

CALI ZEROS (A.b) 

C C A LCl.ll A 1 F PITCH 

c is the overflow fi a& sfU 

tF (JFLAG.F9.0) go to 100 

l THIS SECTION PlcK^ UP TH[ Zt_RO uKOSSiMK FOR AN OVERFLOW. 

„ = 1 ! Y F S , F I K i 1 HE FIRST ZERO CROSSING, 

bn J F ( A ( J ) . Hf . o > GO TO 70 

u = J + 1 

IF ( J.GT.Iphq) GO 10 ? 50 



CO TO bn 
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7H fndi=j ! Found it. 

r / = tIF'pl-Sl ART IC/UCULATt THU PI I Cm PfKIOp. 

SlAPT=LFDl IPlSET THL STrRT POINTER. 

TNDFV=2*J ! CALCHLATi POITJON. 

rUF ( iNOpx ) = <4 U °b ! MARK HIT PITCH FPOcH, 

RU = I NOP X/ 2 r 6 ! T | T S SLClIOM DFT ERr I NCS THE 
] F ( ( KU*?5G) . rO. INCL> ) KU = KU - 1 ! RFC • NO. THAT 

]REr = KU + k\/AR ILCPRLSPONDS TO THE SPEECH FILE. 

IF (IREC.pr.P51) IfUC = 251 

lF (7FLAG.FQ.0) GC f C 80 ! 1 S THIS A S1LFNT AREA? 

T = -I * w lYt'S, CHANGE PI I CM PFRlOp TO A NeG. VALUE. 
7 F L A G = 0 ICIEAP T F ( SILENT FLAG. 

80 |,RJTE< 5* JVAR) h, 1K[.C 

] = J 

v.FLaG = 0 ICLFAR the OVERFLOW FL A G. 

GO TO ?J0 

100 IF ( T . GT . 1?80 ) GO TO 250 !MO, APE WE THRU? 

TH = SUK (IR) ! 0 [ TAIN THRESHOLD. 

IF (TP.fE.U) GO U 120 
/FLAG = 1 
GO TO ?00 

l?n JF ( H ( I ) . L T . T H ) GO T (. 200 !NO, IS I HE PEAK VALUE< TH? 

^=J INC, GO CHECK 7f RO CROSSINGS. 

1 K C )F (A(J).mE.O) GO Tc 180 ! ZERO CROSSING? 

^=J+1 !N(, BUMP THL COUNTfS* 

IF (J.LE.12RO) GO 1C IDO ! THRU BUFFER? 
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START =si ART -1 2tsO lY^S, STORE START LOcATjON. 



I t'r 



iyo 



F LAS = 1 

GO rn ?eo 

[ Ry D 1 =J 

r = tnni -st art 

IN9EX=? * J 



!S(T rut. UVtlULOW FLAG, 

! be RL AD 'HP BUFFER. 

! V L S , ZiRO CROSS J MG OccURRFD, 

! t a I C U L A T t THE PITCH. 

! C / | CULAT t POSITION FOR MARK , 
r U F ( I N IJ F X ) = 409b ! p /■. R K 1 Ht RIlcH EPOCH. 

KU = I N [j E X / ? r; G ! CALCULATE RfcCORU NUMBER 
IT ( (KJ<R5G) ,EQ. Jf’Dl > ) KU = KU - 1 
1 R Ec = KU + kvmR !bf SPEECH WAVEFORM, 

IF (lREC.Gr.P51) IhEf = 251 

JF ( 7 F L /-, G . E 0 , 0 ) GO TO 190 ! IS ThIS A SILENT AREA? 
I- = -l * M ! Y e S » DENUTf WJTH NrS. VALUE. 

7FLAC = 0 

L R I T E ( 5 • J \/ A R ) n, IRlC 

START = FwUi ! R l SET THE START POlMIER. 

I = I + 2n ! R L S E I I H Y MINIMUM PITCH, 

ro TO 230 



2 v r 1 = 1 + 1 

PAP IR = I/25G 

I F ((]R * d 5 A). ME. j) j. R = IR + 1 
(-0 TO loo 

2 bo IF fJFLAG.tG.O) ST A( T = STAR ' - 1 2 €« 0 
0 wRITt Ojt THE ruF^ER OPT C THE FILE 01 * THE DISK, 



r 0 ? 6 0 I A r 1 ,10 

UR1TE(4»KVAR) (RUF( IE f (CIA -T ) *2 56 ) ) * ID=1 i 256 ) 



n rOMTiMiJf 

C rlF-AP THF PJFTEIK, RUF* f C R FHC NEyT CALL, 
r o ?70 ] c=i , pruo 
R U F ( T C ) = U 

2/0 cCNjTJM'Jf 

r o to io 

:v>r rNDFiiE 2 

tNDriLE 3 

fMoftle n 



si op 
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PEKS 




FIGURE A- 6 
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C 

C 

c: title: peks.eor 7 boiio 

c 

c THIS SU 3 POUT jNe; p t M (j S Th[ PLAKS of A w a v r f 0 k M that HaVf. 

C PFEN STOFTP in a RuFEtR t T A 1 B o F i . I HE RET UKN CONTAINS AN 
C ARRAY "F" WItH IHT VALUE pF THE PF AK IF 0 ri E OCCURRED ANU 
C I'A.S G R E. A T r R T H a N PERU* IF NO PtA* OCCURRED aT THAT 

c point, Then a ?ero is returned, 
c 

C FROGRAh PARAMETERS: 

C THE EEERGY SyN'lATJON IN pO-LOOP 3 HAS A SCALE FAC t OR THAT 

c is usfD to prevent overflow. (In ihis case 100. is useu 

L EOR TIE SCALf FACTOR.) 

c 

(• TIC ARRAY "SpMiti CONTAINS THE THRESHOLDS THAT 
c HERE CALEULATEn FpGM THE ENERGY MEASUREMENT. THESE 
C THRESHOLDS ARE RrTUKNLP TO THE. A I N PROGRAM FOP USE IN 

r determining the correct efaks, 

C 

C THE THRESHOLD SETTING TH/\T WAS LINKED 10 a PARTICULAR 
C Mp A S 1 R E v i T N 1 CF ENERGY HAS PETER" 1 INF H RY EXAMINATION OF 
C SEVERAL WAVEFORMS AimP PLOTTING ENERGY VERSUS TRUE 
t non/jL PFMS. THIS ARE; of iHt PROGRAM might HE IM- 
C PROVE F 3Y A MORt D E T A 1 1 F t STUDY OF THE DATA FOR THE 
c OPTIMAL THRESHOLD VAl UES. 
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c: 

C PARflMFTEP "JFI-Ali" JS USEf TO DENOTE A NEGATIVE 
C SLOP[_. TFLAG = 1 - N L 0 A 1 I V t_ SLOPE 
C TFLAG = 0 - POSITIVE SLOPf 

(. 

GU3POUTINF Pf“KS ( E , f , I V AK , SUP ) 

C IMtTALIZA I ION 

IM ret * Fa I Bl'F l ?b6 0 ) , El 1280 ) ,f- < 1 280 ) , SUM( 5) 

COMMON FATBUF 
LO GO I a = 1,5 

C PfPFORFi AN ENERGY SUNFlAljCN OVER °56 POINTS OF THE 

C clottal VAVEFORrt. 

SU V IY = 0.0 

EO ?> Id = 1 ,*>11,2 

SU v IX = SOf'X + ( ( ( pATBUI ((<lA-l>*512)+jB)-2048)/l00)**2) 
? CONTINUE 

C THIS F ART of PROG, ASSIGNS A TRt-SHOLO VALUE BASED ON 
r THE ENERGY ShN^ATTOiV CONflTLD ABOVE. 

IF ( SUMX.NE. o ) GO TO 10 ! SUi“lX = 0 pE A|\|S SILENT AREA 

I HR = 0 
CO TO 55 

IP I F ( SljN X - 12000) 1 5,40,40 

15 tEISUMX - 55^0) 20,35,35 

20 JFIRIJKY - 1 o 0 0 ) 25,20,50 

25 ■ l HR = ??0n 



(0 TO 55 



<0 



lHR = P&on 



3 b 

4 n 
4 b 

bf 

bb 

bn 

C rjMD 

inn 

130 

14 f 
1 b 0 
3 br, 

i r o 



f'O TO b5 

1H3 = 3000 
GO TO 55 

IF ( S U v l X - 1F000) 45, 50 i 50 
1 HR = 3t|0 0 
GO to 55 
T Ha = 3000 
SU'I (14) = JHR 

C 0 M T T WUr 

1H£ PF4KS (-OR THL S L GM L mi of DATA, 
ro ioo kit= 1 i 1200 

F ( KH ) =0 
rOMTIIMUf 
K = 1 

F? = 1 

IF ( FATdUFf nt? ) - FMFUF(K)> 140, 150, ISO 

1FLAG = 1 

GO TO 370 

F = K + 2 

GO TO ljO 

1FLAG = 0 

LO pio h=3,255y,2 

F = 3-2 

IF (FAU'UF(K) .GT.FATrUF(P) I GO TO 200 



IF ( IFLAG.rO. 1 ) GO 1C 205 
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?('0 
2 o *> 
r i o 



JF ( FMBUF( K) ,LT.?Uf ) GO 10 205 

F(*2) = FATBUF* ( K“2 ) iFUUMl.) tHL PfAh 

IFLAG = 1 

GO TO 2 o 6 

1 F L A G = 0 

Y2 - K? + J 

CO'JTTMUf 

R t: T I IR M 

F NO 



ZEROS 




FIGURE A-7 
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C 

C 

c 

c: nil r: Zeros.i-or 750110 

c: 

( this subroutine: ftnos mi zero crossings of a wavl- 

0 FORM SlOPFD jN A RUFFtR. \[ A ZERO CROSSING OCCURS AT 
L SOME FOiMTi THaT P01N1 Is MftRKEU WIlH A ONE IN Tut ARRa t 
r: I'C" . IF NOT, A ZERO is liFT IN THE /\RI<AY. 
r 

c PARAt l-ters used: 

C "IFLaG" IS USED I m this SIHROuTJNE TO ULNOTfc AM 
C A VAu E OF THE DATA THAI jS NEGATIVE. 

C J.E. IFLaG r 1 DATA VALUE IS N t G A T I V E • 

C TFLaG = o data value is positive. 

C 

SUBROUTINE ZEROS ( C , [ ) 

0 

c: initialization 

J NT FGEK CA I B|'F( RbGC ) , C ( 12BU ) , D( 1260 ) 

COMMON catbuf 
ro i n 1 1 1 =i , 1 260 

C ( I T I ) = 0 

if font in J r 

' |> = 1 
In 2 = 



1 



25 I F ( r A T BUM N ) - ?04f) 30,40»50 

30 1FLAG=0 

GO TO On 

C rAUK A ZFRO CROSStN'G. 

4 0 C< M?) = 1 

A = A! + 2 
f,2 = N? + l 

IF (F.BL.2560) 00 Tc 200 
GO TO 25 

50 1 FLAG = 1 

00 IF (CATrUF(N) " 204t> /0»4U,n0 

70 IF ( IFLaG.EQ.O ) GO 10 7 5 

C ( M?) = l 

I FLAG = 0 

75 A = N + 2 

A 2 = A'? + 1 

• IF (F. Of. 2560) GO T c 200 
GO TO bo 

00 IF ( IFLAG.LQ.l > GO TO 75 

0 ( M? ) = 1 
IF LAG = 1 
(0 TO 75 



20 0 



RE 1 URN 



- 171 - 



ADDNS 




FIGUR E A- 8 
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C 

C 

C TIT|.F: ADPNR.FOK 7501P0 

c 

C TUTS |-R0r,R/\' v l ADDS Two FIlFS T ogethfk WITHOUT a SCALL 
< FACTOR. IT is .v'E/’Nr TO Lt USED ON An ODD - F.VEn FORMAT. 

C (I.L. 0\Jr F J t F OQO FOKFat , THE OTHER FVEN). IF SWIT C H(U) 
C IS UF T Hr PROGRAM ALLOWS THE INPUT aF'D -OUTPUT FILFS TO 

c pf n a ffs. jf switch (0) is down, thfm thf program 

C ASSURES the DEFAULT NA Fit S that Apt. USED IN SERIFS WITH 
C THE CCFPl FTt SHI El SYSifm. 

C 

C SUPRCLTlripS USED 

C I SSW I ( SEF REVArp FOR FX[ [ ANANT ION. ) 

c 

c INITIALIZA I ION 

INTFGFR DA | Bl'F ( 2bG0 ) , CATBUP <2560 ) . BAt?UF ( 2560 ) 

TATA IVAR. JVAR«KVAK/?,2.2/ 

LRITFI 7,1 ) 

1 FORMAT (« • , * PROGRAM CONTROLS - - - SwUCH (0)') 

URJTF <7,?) ' 

2 FOR^aTI • •» ’SWITCH! c) DOWN - DEFAULT FILES USED*) 

LRITF <7*3) 

3 r C R m a T ( . . JNPUT - SIGNAL.SSS*) 

' WRITE! 7, q ) 

R FORMAT!* * * • rU T PUT - SIGNAL. PAT*) 
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UklTfM 7,5) 

B T(« • • « c W I T C H ( q ) UP - I/O FILLS MAY BE CHOSEN’) 

l AJsr ’PAUSE* 

IF ( ISSW'I ( 0 ) .LU. 1 ) (. C TO 6 
CALL ASSIGN ( 2 , * KK 1 • S I GNaL . SSS ’ ) 

CALI ASSIGN (3**kKl;PIICH.UAT*) 

CAu ASSIGN ( 4 « ' K K 1 : S I G N A L • 0 A I * ) 
f-0 TO 10 

6 WRITE (7,7) 

7 r 0 R Pi A T ( • » * * I* A ME OF INPUT HLF <F I LN A« . EXT > » ) 

WRITE ( 7, A ) 

A FORMAT <’ * »• ') 

C ALI ASSIGN ( 2 » * R K 1 ; SIGNAL. SSS* ,-14) 

CALI ASSIGN (3. ’KKi;PITCH.DAT* ) 

WRITE (7.9) 

«•) FORMAT (• ’.’f'AUE OF CUEPUT FILL <F 1LNAN. ['XT >’ ) 

WRITE( 7, A) 

(ALl. ASSIGN (4. • KK 1 : S I GN AL . D A T ’ ,-14) 

](J f'EFlML F 1 1 r ? (2tJl.c5S»U»IVAP) 

CLFINL file * (2t>l.j:S6.U» JVAP) 
rEFjNE FILE 4 ( ?U] , 256 , U , KVAH ) 

( PROCESSING BEGINS LlLRL. 

11 jF ( IVAR.EQ. 2^2 ) GO if 30 0 

c fill pUFrrRs. 

CO 20 1 A = 1 . 1 0 

REAP ( 2 ’ IVAR ) (OA)PUf(JB + ((TA-1) * 256 ) ) , IR=1 . 256 ) 
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20 CONTINUE 

r o 30 i c = l , l r 

I EAP ( ? ’ JV AK ) (CAlRUf(IO + ( ( I C - 3 ) * ? 56 ) ) » ID = 1 » 256 > 

30 L C NlT I A'Uf 

C APO IRE. Tl'C TOGfTHEH. 

TO 40 J E = 1 , 255U.2 

RATpUF ( IE ) SDATWUK ( H ) + E A I B 1 IF (IE) 

F-ATpUF ( IE + 1 ) = CATBlF'(IE + l> 

40 COMTINUL 

( WRITE THF OUTPUT 

10 30 IF = 1,10 

WRITE (4'KV/AP) (UATlUF(IG + < < I F - 1 > * 256 ) ) , I G = 1 » ?06 ) 
30 COMTTNUf 

GO TO ll 

300 f NOplLE 2 ! Tl 'RW I H 1 f G S OFF. 

ENDFILE 3 
FNDFTLE 4 
STOP 
FNO 
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REVAMP 




FIGURE A-9 
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C 

c ini r: Revamp. fui-* 75U1.30 

c 

C SUBROUTINES USED: 

(' TSSW1 - THIS SUBROU1INF 1? WRIT I EH IN ASSEMBLY 
C t ANGUAGE A|\|U Api OWS FUR HE DISPLAY REGISTER SWITCHES TU 

c rt used as program controls. 

c 

L r OR su ITCH ( 1 ) UP 

C T R PUT AMD OUlPlJI EIlES M/ Y BE SPECIFIED. 

C 

c for si- itch (o) noi.'iv 

C INPUT - SIGMAL.IJAS 

C OUTPUTS - 1. SIGNAL. DAT 2. P11CH.TAU. 

C 

r THIS p K 0 f R A V) IS used 1C REVAMP A SIGNAL FILE THAT HAS 

C EFEtv edited, the OuP v/ l ues of the signal are preserved 

c IN THIS TPA'MSFORvaTION. l HE EVEM VALUES OR HITCH NARKS ARE 
C 7EROES JNLESS THEY APE Cl TO 40/5. AT THE SAME TIME 
C thf: hi stance between pitch marks IS MEASURED and 
C STORES IN PITCH. TAB FOR USAGE I N OTHER PROGRAMS. 

C 

C INITIALIZATION 

INTEGER DA I'BliF ( ?bGl! ) ♦ BATBUM 2560 ) « STAR 1 « A ( 256 ) 

[ATA SlART»A/2«?5f * 0/ 



LATA JVAR, JVAR,KVAR/?,2,2/ 
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C. PFVlfU PPFVIQUS PTTCH. T A[. FILL 

CALI ASSIGN (?. »HH1 :( I ICH. I AR* ) 

[EFlNE FILL P (1 UOU , ? . U» JVAR ) 

1 HEAP (?• J\/AR,Elvin = L ) I K EC 

IF (P.LT.O) MlRLC) =1 
CO TO I 

? ENDFILE 2 

C NOW RF 1 MIT I Al.l7t . 

IaRJTF (7.3) 

3 F ORf’AT (’ • , • SWT 1 CF ( 0 ) DOWN - DEFAULT FILES USEp. ' > 

WRITE (7.4) 

4 p 0 -I m A T ( « » ,MNPUT - SIGNAL. PAS; OUTPUT - S I GMAL # DA I * ) 
W R 1 T r (7.5) 

5 FORMAT (’ • , * SW I I CH i 0 ) UP - I/O FILES NAY BE CHoSElM') 

I- A U S F ’PAUSE' 

IF ( ISSw. I ( 0 ) .EO.l > l C TO 6 
CALI. ASSIGN (2. * KK1 ; SIGNAL . DAS* ) 

CALL ASSIGN (3* 'KK1 ;SIGNAL.DAT' ) 

CALL ASSIGN (4.'HK1 JPITCH. IAR» ) 

(0 TO 10 

0 WRITE (7.7) 

7 FORMAT (• ' , * N ANL 0( INPUT FTLL <F I LNAN . EXT > ' ) 

WRITE ( 7 . P ) 

A |- 0 -0' r A T (• •) 

CAU ASSIGN (?. 'KK3 . ; SIGNAL. DAS' ,-14) 



WRITE (7,9) 
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9 FORMAT (* •,»NARt G[ OUTPUI FILE <PIL.NAM.EXT>*) 

V R 1 T F ( 7 * 0 ) 

CALL ASSIGN (3» »KK1 JSIGNAL.DAT* 

CALI ASSIGN (4» M<K1 ;PITCH. I AH* ) 

10 I E FJNF FILL 2 { 2b3 , 2 5 6 , U , I V AR ) 

CEFIME file: ? (281 ,256, J, JVAR) 

FEFtNL FILE a (100C,?,UiKVAK) 

„VAR = 2 

c are irE thru? 

11 IF ( IVAR.so.252) GO 10 300 

C FILL I HE PUFFER WTTH UA 1 / FROM INPUT. 

L 0 ? 0 1 = 1,10 

EEAP( ?' JVAR) ( DAT BUf ( J + ( < I - 1 ) *286 ) ) * J=1 • 256 ) 

2 0 CONITINUf 

C TRANSFER THE odd c aEiRE-ES as they are. TRANSFER the even 
c SAHEL F S AS EERO UNLESS T ) F Y EQUAL 4095. ALSO THE 
C riTCH.TAp F 1 1 E is BEING LRITTEN OUT IN THIS LOOP. 

LO 50 U = 1,2859,2 
PAT OL'E ( II) = DATBL'F ( J1 ) 

IF (TAIBUFI 11+1) .CL. 4075) GO TO 25 
P A T P U F < Jl + 1 ) = 0 

CO TO 6C 

25 PATPUF ( Jl + 1 ) = 4095 

E = ((Il+l) - START)/? 

,.U = ( U + 1 )/256 

IF ( ( JU*?5S) . rQ. ( 11 +l ) ) JU = JU -1 
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uKt-C = JU + JV/«R ! (. / l C U L A T t REC. NO. 

C CHECK FOR S1LEHI cRLAS 

If- (R.LT ,200) CO TO <(0 
IF ( A( JPEC) .HE. 1 ) Cc TO 40 

K = R * -1 ! AREA IS A SJLLNT ONE. 

40 In RITE ( 4 ' K V A R ) R» URFC 

START =11+1 
bo COVTIHOe 

r. K RITE THE OUTPUT INTO SIcfAL.UAI FILE 
EO 70 1? = 1,10 

RITE ( 3* JyAR) (RATB0F(I3 + <( 12-1 ) *25&> )« 13 = 1 . 256) 
70 fOMTJMUf 

C PFSET ST ART 

start = Start - 25 h 
GO TO 13 

30p ERDFILE 2 ! TURN THlfv.CS OFF. 

ENDFILE 3 
ENDFILE 4 
STOP 
E NO 
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I TILL ; SHIELD. OOc 
oocumentation fok s h v o 3 add shv04 



APSTRAC T 



T 1 1 E S L F’KCGRf v l c ' Il u iPL| RFHT IMF SYSTEMS DESIGNED AND 
MTeLAIFD PY V A UE ri SHILirS IN AN MIT MASTER'S THESIS FROM 
1X70. THE PURPOSE OF THESl SYSTEMS WAS TO S:PARATF A SIGNAL 
FROM NO 1 S E (USUALLY ANOTHl.lv SPEa*lR). THLRE ARE TWO SYSTEMS 
IMPLEMENT EP, AND 1HFSE DlFlfR IN THE MANNER OF UNVOICED 
SPEECH PROfTSSlNG. S H V 0 7 USES I HP ATTENUATED input METHOD 

WI'TI F SFVDU USES IHF i NER T j A L F I L 1 ER 1 NG METHOD. 

SUBROUTINES USfD 



i. window ( ntype,k) - this subroutine allows for ihe 

/ 

USrP To PRESCRIBE THE WINLCP FUNC'IPN DESIRED IN THE PROGRAM 
Arp TO rEbIGMATE IMF VALUE pF K FUR THE FILTER LENGTH. 



2. CnFF (NTYPT ,K, A.L ) - This SUBROUTINE CALCULATES IHE 
LOrEF I c I t MTS Of TUT FUTlR PRESCRIBED IN THE WINDOW 
SUPROUl THE AND SjORrS THESl VALUES IN AN ARRAY "A" THAT GAN 
HI OF MAXIMUM DIMENSION 01 IS. 

0 1 pfR f ROGRaNS USlq JIM CONJUNCT 1 0 WITH SHIELD. 
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l. ^ r s r w i p - nils program sa »ipi.e s the analog wavefurm 
A 1 1<J KHZ. OK Each OF ? CHANNELS (4 aNU 14) AND PLACES ITS 
OUTrUT c l\j FILE SlbNAL.VXX . the chammels are STORED in 
J HTERI-F A V 1 HG FORMAT, 

C H A N M r L (4) “ ODD SAh| | E VALUES. (ODD FORMA!) 

ci ./ifi \ ti. (i0) - eVtm Sample values, (even format) 



2. Sf Pr< - I HIS PROOF A M SePaRa ' fc S THE TWO CHANNELS ON 
0 N F FILF AND p I ACFS EACH G( A SePaRaK FILE. 

inf l i - signal. xxx. 

OLITPJTS - CHAMfEL ( A ) - SIGNAL. SSS 
CHANNEL (14) ~ PITCH. SSS 



BOTH APE ASSIGNED 1 HE FORM A T OF A CHANNEL (4) ONLY 
Mir. ( I , F , ODD SaMIPLlS - DATA? EVEN SAMPLES » 0 BY 

m fault.) tf oeStrep hie Second output can re put into even 
E or p AT . 



5. 


^IT VO? 


THIS 


| POGRAM DETFC1S 


pitch 


FROw| 1 HE 


PITCH. SSS 


F I Lf . 


THERE 


/ pc 


TWO OUTPUT 


FILES 


FROM THIS 


PPOCRAE . 


FTRST , 


PI TCH. • At IS 


a tabular file 


OF 


THE VALUES 


OF THl 


p TTCH 


PLRTOu. 


T) IS 


IS USED AS an 


INPUT 


to shield. 



SECOND, PITCH. PAT is A Pill ihat ts used 10 MARK PiICH 
PlFlOUS in A WAVEFORM F I L | , PITCH. PAT H AS TWO VALUES ONLY. 
(U OR 4095) , THE VALUE 41<E IS USED TO DEMOTE AN ENpPOlNT 
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oi- a pitch prRior. this is ijssp as an input to adohs. 

PITCH. DAT is I\j EVEN FORMAT. 

'+. Annns - THIS PR( CRAM ADOS T O f ’ET HE R TWO FlLt.Si 
SJr-f'AL.yXX A IM u PlTcM.UA I . SINCE UNt E J L L * SIGNAL. XXX, IS IN 
Opr. POP/, T /»N0 PITCH. DAT IS> IN EVEN, I Hf CHANNELS are kept 
SEPARATE but ^UT i\)TO ONI (TIE CALL TP SIGNAL. UaT, 

INPUTS TO SHIELD 



1 . 


SI G N A l .UAT - FROM AOONS 


2 . 


FITCH. TAR - fROM pH VO? 


OUTPUTS PROP 


shield 


1 . 


OUTPUT. ua i 


At GORITFP 





1. Tuf 


PROGRAM FTRST PROVIDES THE USER PITH THE CHOICE. 



OF TYPt OF WlMUOW AND THE f A R AML ILK K. ( 2K + 1 COEFFICIENTS 
APE USE! IM TH[ FILTER)* 

2. I HE RROrKa 1 ”' THEN caLCULAIES THE cOtFFICIENTS FOR 



T Hf FILTER and STURTS THEM JN ARRf'Y "A" 
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4. The NeXT Srcliun / I LOWS ROP TrE inii ialization Op 
1 HF LOGICAL OMjTr TI’aT will RL USLD IN This implementation, 
a. LOGICAL Uf;TT ? - LSFU FOK THE SIGNAL FILE. 

D. LOGICAL U!\JTT 5 - L? EU FOK THt PITCH FILE. 

C. LOGICAL UNTT 4 - USED FOK THL OUTPUT FILL. 



NOTE: the INPUT aNU OUTPUT signal FILES, on LOGICAL 
III' 1 ITS ? aF’E) 4 AKE D 1 V 1 1 > t L INTO RECORDS OF LENGTH 256. I HE 
R ] T C El FiLt is P I V I D F R TNTO RECORDS OF LENGTH 1WO. 



4 . 1 H E BUFFER Is THEN S LT Up INITIALLY UY LEAVING IK 

records full or zeros <ir is computed in the following 

MaMVFR ; ( 1/9 Op I'HF WINLCW LENGTH) IN SAMPLES * 1/256 

(RT COKL/SA^pLEs) + 1 I F f THE NUMBER OF SAMPLES IS NO I A 

M'l TIPLf Or THEN THt SIGNAL IS TaKEN FROM T HF FILL AND 

pi aced into the puffer. i he numpfr of records read is 

VAPIABLF OrPFNpING ON THF fJRST VALUE OF PITCH AND ON «. 

5. IN THIS PROGRAM I)F VALUE I IS USED AS A BUFFER 
POINTER. THL VAlUF 1 fOlNTb '0 THE CENTER OF THE FlL'CK 
THAT IS BEING USED. I IS COW SET UP TO BE THE FIRST SAMPLE 
IN THE FAT A SECTION OF I HE F UFFER • 

6. IMF V A LllL OF riTC| IS READ AMD THROUGH THE SAME 
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SI: PITS OF MtLASUKErrFuTS '«S FOUwO In THF I N I T I AL 1 7 . A I T I OIM 
FIT CESS, 1|IF ALGORITHM Cl'UKS TO Sf E IF 1 «E t XTRt TE VALUpS 
Fl-FPEIJ IN THE f I F K F P F Im F L LGUaTION APF within the buffer. 



7. IF TH[ Y aRF HOT, 1 HE BUFFER IS SHIFTED; THE VALUE 
OF I K[" p OS 7 T I Of En A f 'D F't-l | A,TA IS READ IN FORM THE FILf. 



8. MOW The CALCULATIONS cam be perfomed, i he 
U lFFERtfCE EQUATION I? J N{ L E^E|\TEU A r 1 B An oUIPUT POINT Y IS 
CA ICULATEJ AMU WRITTEN ONTc THE OUTPUT FILE. 



9, T|,iF PROCESS CONTll-UFS UNlIL > h'E HITCH FLAr IS SET 
BY PLACING TPE Ma«K in THE IMPUT FILF. A NEW VaLUf OF PIlCH 

period is them read, aN'u tif filter is modified for this 
change. 



10. A SLIDE OPEkATIOi. OCCURS WHEN THE VALUE OF I IS 

iuo stall to accommodate hf current filter, the buffer is 

KE POSITIONED AND 1 RESET Sc THAT ALL INPUT VALUES NEEDED BY 
IHF FILTER ARE P R E‘ S F N f l Y I(. THE BUFFER. 



11. THE PROCESSING C G f T I N U E S t 'NT I L A SILENT A«FA IS 
OtTFRKlNFO, A\Iq tMEN Af I E K FILING THF OUTPUT FILE WITH ZFKOS 



IHF PROGRAM FIAlTs 



SHV03 




FIGURE A-10 
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SHV03 
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SHV03 




FIGURE A -1 0 (CON'T.) 
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SHV03 



2 70 




300 




End 



400 , 



Quit 



FIGURE A-10 (CON'T.) 
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C 

C 

C 

C 

C 

C 

c 

C 

C 

C 

C 

C 

C 

c 

c 

r 

c 

c 

C 

c 

c 

L 

c: 

r 

e 

c: 



iinr: r h v o 3 author: r. fk/vilr 

ADAPTED F h TliESJS PY V. SHIELDS 

this [ F\ o n r a 10 uses thf pm pod of shiuds to hep for P'1 

THE SF F Ap A T I ON OF A SPEYER Af g D NOISE UY r,EAnS 
OF A DIGITAL C 0 hi B Fll TE.K. 1HIS PROGRaP'1 USES THE 
ATTLFjIATTON KETHOn FOR U| VOICED SECTIONS. 

SURRCLUrTS USED: 

hi mode ( ht y re » k ) 

COEF < n t ype ♦ K * A » l ) 
fa.rapttlps used: 

ARRAYS ; 

1. OAIbUF _ contains the INPUT SAPIPLE values, 

R. A - CONTAINS THl COEFFICIENTS. 

3. OUlpUT - TEMPORARY BUFFER FOR OUTPUT VALUES. 
va r iafle9; 

1 . SFLaG - USED TO INITIALIZE OUTPUT FILE 
SO THAT STEP DISCONTINUITY DOES NOT 
OCCUR. 

? . ZFLaG - OSLO TG DENOT HFG. OF SENTENCE. 

3. PFLaG -USED TO | t NOTE Np-' PITCH PERIOD. 

< 4 . LARGE - USED TO PENOTE UNVOICED SECTIONS. 

F. COM - 1 Hr ATTFfLATlON CONSTANT FOR UNVOICED 
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i' sections. 

r 

C INITIALIZATION 

c 

Ll v IFNSloN M'5) 

1 N 1 rHE « DA I Bl 'F c 4096) , START » OUTPUT ( 256 ) .ZFLAG.SFLAG 
logical PflAG.LAKGE 

[ ALA DAyCut- , c-flAGi Ol TPUT/ tf 0 96* 0 , 1 1.256+0/ 

LAI A IVAR,vJVAR»KVAk/2i2,2/ 

C CALL THE FOLLOWING SUBROUTINES POW I N I I I AL I Z AT ION . 

LR1TF( 7,1) 

3 FORMAT ( ' '.M/ALUL Of A I TEN « CONSTANT? <1»< F3.2>') 

aCCFPT 2, CON 
2 F 0 R N A T ( F 3 . 2 ) 

CALI ASSlGN(? t 'RKlIijlGNAL.lJAT* ) jOlJTPuT FROM ADONS,. 
CALI ASSIGN! ‘PITCH. TAB* . 9) IOUTPUT FKON PITDLT. 
CALI ASSIGN(«» •OtnPLT' «6) 'OUTPUT OF THL COW 

! FILTER. 

rEFTME. Fllf 2 (253 i 2 5 6 * U , K V A R ) 
r EF INF file: * (1 000,2«U| 1VAR) 

C WIMOCL OFTt RFl.MAT 1 ON. 

CALI W If, DOW <MYPF,h) 

C CALClLAIF coee ficte nts. 

call CUFF ( NTYPF » K 1 a » L ) 

C SPECIFY the LENGTH OF TH[ outpui fhf. 

I I "ITT = 253 + ( K * j ) +1 
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ILFTIVE FILE II (Llflll .26S»U« JVAK) 

( FILL THE PUFFER INITIALLY 

C RFAU I HE FIRST V A (Ut FROf THE PITCH TABLE FUR 

c inn ali?ai ion Supposes. 

20 L L An ( 3 • I v A R ) M 

c rOUBIL 1 H I s VALUE TU CTMf FNSA1 E FOR TWO CHAU. 
[ - = I A o S ( P * V, ) 

c is n too larce? 

1 P < *•". Gf_. 40 0 > 00 TO ?0 

c row wt V.y set up duff i r , take "ieasurei'tlnts. 
be c = L - 1 

i v = ( ( c / ? ) * m + 1 

lR = J\//2[ife 

IF < JR*?5A»NF.1V) II =IR + 1 

•j F-TKi ?56* 2 

I=( IP + 1) - (K * M) 

0 FILL THE PUFFER FPOli WHERE 1 HE POINlFP STARTS 

to F D 1I = (2*TR)+I«lt 

RLAP(P'KVAK) (UA1BU|(J + ((II-1) *256)), 

on comtinue 

I V A R = ? !K[SEI PITCH. TAB FILE. 

C row of Glr THf P«OCE sSir'G. 

1 op OF lac- = .FALSE. ! T » . 1 S FLAG DPI. WHEN /E 

?f LAC = 0 
R E A n ( 5 • I \/ A R ) H 



J=1 , 256) 



CHANGE PITCH. 



J F ( r. LI . p ) 7 F L A G = 1 
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f = i r n s c ? * v, > 

C T S M TOO l M c E 9 

IF (M.LT.400) GO TO ] 30 iM>400 COKHES. TO l J ITCH<50. 

C YT S, 1 U r? M OFp p 1 1 TEK. 

1?0 LARGE = • T K U F • 

T I F ( PFL T f, ) bo TO 10 0 !IS PE I- A G SET? 

c pflag is noi srr. 
r APE WE 1HPOUGH with thf e offer? 

IF (I.GT.409F) GO Tr 127 
< MO, GC RfTUKm AM OUTPUT ( PINT 
ASS TON ]23 TO KXX 
f 0 TO 2 0 0 

0 YE S » W E ARE THRU PUFFER , READ JN (I IMEW BUFFER , 

1 ?7 |0 )J>6 JJ : 1,16 

IF < KVAR.GE* 2 C 2 > GO 1 0 300 

RE Ap ( 2 • KVAR ) (DATPUfdCB + ((JU-1) * 256 ) ) , I CB=1 , 236 ) 
1 2 P CONTINUE 

129 i = i 

GO TO 125 

t M WAS NOT TOO L A R P T » PROCESS THE SIGNAL. 

130 LARGE = . F o L c E . 

G TAKE fEASUREmEMTS. 

JV=((C/?) * ME + 1 ! IV = l/2( LENGTH OF T HE WINDUW) 

lR=TV/2f G 

IF ( IE*?56.i\iF. IV) lf« = IR + 1 



1P = IP* 25fc 
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C IS T I R P T T C H FLAG SlT? 

1“C Ih(PFLAG) bO TO 100 

c r-o, rcos the: puffer heed to be shiftfd? 

1 bo J R ((I + lV) t LT.ROVt) GO TO 160 

C YES* SHJrT THE RUrFER 

3 H 0 1 n=I 

KVARrK 1 /aR- ( ? * IK) 

[0 1 5P JJ=1 , 1 6 

IF (KVAR.Gt.P52) GO TO 300 'ARE WE THRU? 

REAP ( ? * K V A R ) ( OAT E'LF ( J + < ( J J- 1 ) *256 ) ) , J = 1 , 256 ) 

1 5r C OY'TT M’Jl 

c reset thf POINTER i. 

150 START = (lb -(2 * 3 1 ) ) * 236 

I = I HOLD - START 
C-0 TO 1 90 
r is I TOD SMALL? 

16n IF(T.GE.IV) GO TO 1^0 
(' YESi SLIOE. 1 (IE BUFFER. 

170 ^HOLO = I 

^VEx = IV - T 
^VEy = JVFX/P56 

]F((JVLY * ?*6) . HE .^VEX) JVEY = JVEY + 1 
KVAR = (KVAR - 16) - JVEY 

1.0 175 IDy = 3*16 

REAP(?’KVAR) (uATFUf ( IDX + ( ( I pY-1 ) *?b6 ) ) » IDX = 1 t ?t>6 ) 



175 



( OMTINUE 



- 194 - 



c: FT S E 1 I 

T = JHUlF) + ( 2b6 * j\/EY) 

3 9f ASSl<TJ 140 TO ft XX 

( roMPin ai tom /\nd output Section 
pop 1ND = I 

? To IF ( TMP.LC.P5f ) 00 lc 220 

1I\1Q = I (V D - P5b 
GO TO 2] 0 

PPG T F (LAKGF) Go 10 ?b o ! LARGE SET? 

0 NO, CpLCUlAlf AMD 0UTPU1 POINT. 

Y = 0 . 0 

TO 24 C NMrl.l - 1 
1 P 0 1 N T = 1 - ( (K + 3 - MM) * N) 

Y = r + ( A ( M I4 ) * PATbLF( IPOI'MT ) ) 

P'+f CONTINUE 

c I? TF IS THE FIRST OUTPUT TO BE WKITIEN. 

IF' (SFLaG.EQ.O) GO if 2 55 
G YES, c c N 1 T ALLOW pRASTIC STEP INPUT. 

EC 254 IU r 1 , ( INF -j. ) , 2 
OUTPUT ( IU) = 2040 
OUTPUT ( IU +■ 1 ) 0 

PH4 CONTINUE 

SFLAO = 0 ! F n R L V L P . 

0 VP I T [. 1 Hr OUTPUT POINT. 

PPG • GUI PUT < INQ) = Y 



(-0 TO 2 7 0 
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2f 0 lF(7FU'G.rO.N GO 1( 265 

CUTPUT ( 1ND) = c04« 

CO TO 270 

?fcr OUTPUT! ]NO) = (CUN * (UfllBUF(I) - 2048)) + 2048 
2 70 CUTPUT C J h* n + 1) = 0 
c ARE Ir F READY To Olir'IH OUlll'T? 

1M TOL.M. ,255) GO H 280 
C YES.UE APE! 

UR ITT ( 4 « JVAR, ElMD = 4 0b ) ( OUTPUT ( KKK ) , KKK = 1 , 256 ) 

20 0 1=1+1 

IF< PATrfuri i ) .F0.4065) PFLAb=3 ‘CHANGE PITCH? 

1 = 1 + 1 

29n CO TO *XX, (328.180) 

C THIS .SECTION IS (jSFU I P fJLL LA • BUF WITH ZEROS WHEN 
C CALLEL. 

300 LO 310 lLAX = JJil6 JDtTERMNL HO*' MANY 

EO 305 JLAX=1 ,256 

E A 1 P U F ( J L A X + ( ( 11 A> -1 ) +256) ) rp ! RECORDS NEED 

in PE 

30p CONTINUE IAUDFU HUH ZEROS, 

31 n CONTINUE 

IF ( I ARSE. ) GO TU 129 
CO TO 159 
C YES, c U I T 

4 0 0 • f N0E1LE 2 



f NDFlLc. 3 



ENDriLE t) 



STOP 



c 
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C 

lllir: SHv /04 AUTHOR: R. FRA7ILR 

r ADAPiru From THESIS E Y V. s^lElDS 

C 

c TI TS I ROoPA"' USES THE METIOu of SHlTEpS TO PERFORM 
(' THE SEPARATION OF A SPt f y \ R AjVjO MO IFF BY i^EAnS 
0 OF A t; I G T T A L COM 3 FILTER. I MIS PROGRAM USES THE 
C ATTEfvlAlTOM METHOp TOR Uf> VOICED SECTIONS, 
c: 

c suproltimfs usfo: 

c HINOOU (MTYPE«K) 

c: ooef (NType»k«a«l) 

c: 

c parameters USED: 



c arrays: 



c 

c 

c 

c 

c 

c 

c 

c: 



c 



c 



1. DA I BUF - C 0 N 1 A 1 [ S THE INPUT SAMPLE. VALUES. 

?. A - contains the coefficients. 

3. OUTPUT - TEMPORARY BuFFTR F OR OUTPUT VALUES, 
VARIAF-LES: 

1. SFL/>& - used to initialise output file 

SO THAI STlP DISCONTINUITY DOES not 
OCCUR. 

?. ZFLAG - IISLD TO f F N 0 T E BF. OF SENTENCE. 

• 3 . PFLaG - USLD TO rENTOE NEW PIICH D FPIO(j. 

A. LARGf - USED TO CENOTE UNVOICED SECTIONS. 
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C 

t JNITIALIZAIION 

c 

l.lflrNSloN r ( i 5 ) 

1NTFGEK DA I B"F ( 4U96) , S T ART » OUT PUT ( 256 ) ♦ ZFLAG.SP'LAG 
I OGtcAL Pfl A n , L. AHCL 

LAI A DA T b u F , cFlAg.UlTPUT/ 4 096*0,1,2564(1/ 

DATA JtfAR * JV AK » KVAK/? » 2 * 2/ 

C CALL THE F 0 L L 0 W 1 N r SUBROUTINES FOR I N I T I AL. I Z AT I ON . 

CALI ASRIGN( o, *PK1 rSIGNAL. DAT’ ) ! OuTPuT FROM ADDNS, 
CALI A Ss I G i\! ( 7 , ' P 1 T L l , TAB’ « 9 ) ! OUTPUT FROM P1TV02. 

CALL ASSIGN(U, •OUTPUT* » 6) ! OUTPUT OF THE COMP 

! FILTER. 

rCFlNE. Fll E 2 (2bl,25G,U,KVAR) 
rEFTME FILE' * < 1000 , ? , U , I V AR ) 

C WINDOW OFT E Rf. I M A T I ON • 

call window ( f’T r Pt , K) 

c CALCliLAIF COEFFICIENTS. 

CALI cot F ( rJTYPL 1 K « A , | ) 

( SPECIFY THE. LENGTH UF TH L OUTPUI F I L f . 

L I "ITT = 2 51 + (K 4 % ) + 1 

rEFTHE FILE' n (I 1M m?56,U* J'/AR) 

C FJLL THE PUFfFR if.'j IIALLY 

C PEAIJ THE FIRST VALUE H’Of THE PITCH TABLE FOR 
C INT TIALI 7AL JON Pu^PO^LS, 

2 0 



F F A p ( 3 * I v A R ) LT 
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c rowBif this \/ai<-'e to u r i i nsai e fob too channfls. 

K = I A R S ( ?*''!) 

NHOLD = N lyl'IS WILL RE USED FOR uNVOICFD, 

C TS M TOO | /iK ( ,F? 

IF (N. OF. 400 ) bO TO ?n 

0 NOW Wt v l a Y S t T up pUFF L R , TaKE ’TEASUPfMLNTS , 
bO C - L - 1 

iv = ( ( c/?. ) + r) + i 
lH=JV/?b6 

IF ( TR*25F,.NF. I'M 1 1- = I R + 1 
jP = tr*? 56 *?' 

} = ( T P *■ 1) - ( K * N ) 

0 TILL the pUFF f R FROM WFLI f THE ROINTpR STARTS, 
ro PO I I = ( ? * T R ) + 1 . 1 c 

FE‘ftr»( ? * KVAK) (UATRUMJ + (tll“1) *?5G ) ) » J=1 , 25G ) 
OP COMTINUE 

)VAr = 2 ! Rfc SET PIICH.TaB FILE. 

C 

C NOW PtPlM THE PR0PESS1NC7. 

4 

1(10 1 -FlAF = .FALSE. ! 1 F J S FlA k DEI. WHE|\) WE CHANGE 
I'TTCH. 

7 F L A G = 0 
FHO| P - N 
|> E Fr> ( 3* I VAR) N 
‘ ]F ( N. LT . 0 ) 7 F LAG = 1 



F = J A R S ( p * v) ) 
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C TS M TOO LARGE? 

IF ( '".Gfl.ttno) m - ft I old 

L TAKE FEASUREFFNITS. 

] V = ( ( C / 2 ) * P> + 3 ! JV = i/2(|.ENGTH OF THE WtNOUW) 

IR = t\//?56 

IF ( I R *" ?5 F • MF , i V ) I (\ = I R + 1 
ll’ = IR* ?5 6 

C TS T|ir P TTCH FLAG slt? 

3. '4 0 IF (PFLAf 7 ) GO TO 3 00 

c: 4 0, [jCLS THE PUF'FFR ['LED TO BE SHTFlFD? 

3 bn H ((1 + IV/),U,'4 0^t>) GO TO 1GO 
c: yfs» shift «he pufflR 
I HOI. 0= I 

KVA<v = KV/AR- ( ? * IK ) 

FO 150 JJ=J , 1 f, 

IF (KVAR.Gl.?52) GO TO 300 ! ARE WE THRU? 

REAP ( ? « K \/ A R ) ( OAT bL F ( J + < ( JJ-1 > *256> > » J=1 » 256) 

150 (OMTINOF 

C reset I 

35° ST APT= < 3 6- ( 2* IR ) > 4 ?5 6 

} = iMnL.o - start 

GO TO ?00 
( TS 1 TOO SHALL? 

U 0 IF( T .GE, JV) ro TO 2( 0 

G YFS, sLIoF IRE RUFFlP. 

3 / n oHOLP = I 
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v. V t y = j v - t 

= jVFX/256 

lF<(jVtY t ? «= 6 ) .Nf . C VEX) JV/£Y = JVF.Y + 1 
kVAR = (KVAR - 36) - JVEY 

f'O 175 J D y = 1 i 3 6 

RLAP(?’kVAR) (UMPUj ( 1 D>; + ( (irY-l)*?56) ) »lDX = l«2t>6) 

17? f 0 MT I MUt 

c reset i 

I = JIHOLD + (256 * cVE Y) 

C rOMPUTAlTON APd 0"TPU1 S t ( 1 ION 
pro IMO = 1 

230 IF ( TND.LE .25* ) GO 1c 220 

I ND = 1 N 0 - 25o 
GO TO 2 3 0 

C PO, CALCULAlt; AMD OUTPUT POINT. 

220 IF ( ?FLaG. ECO. 1 ) GO l 0 253 
Y = 0 .0 

r o pm o pr=i , i -l 
1 P 0 1 N T = I - ((K + 1 - MM) * M ) 

Y=r + (A (MM) * ( UAT LPF ( IPOINt ) - 2048)) 

240 fOMTlMUf 

C IS TpjS THE r IRST 0UTPU1 TO BE! ^RTTTfN? 

26? IF (SFLAP.L'Q.O ) GO -| 0 2 55 

0 YESi | OV T a | LOW PR AS 1 I C STEP TRANSIENT. 
r,0 ^54 iu = 1 , ( INI -2 ) % 2 



OUTPUT (TU) 



2048 
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CU1f>UT (T.U + 1) = 0 
P 0 (| CONTINUE 

SFLaG = 0 ! FOKLVLP 

c i-'RH'l l Hr outpjt point, 

255 IF < ?rLAr . FP. 0 ) 0(jTF L T(IND) = IF I X ( Y + .5) + 204fl 
JF (ZFlaG.Fq.I) OUl | l)T (INU) = 2048 
27o JF (nAlfjUFd + 1 ) .bL .407b) OMTPlITdNp + 1) = 4095 
JF ( P A I R U F ( I + ]>« I 1.4075) OUTPUT ( I Mp + D = 0 
C /'Pfc. Wf 4FADT TO OUPP opt | UT? 

TF ( TNP.pt .255) GO u, 280 
C YFSiWL APT! 

UR1TE< 4» JVAR,ENP = At>0 ) ( OUTPUT (KKK) .KKK = 1,256) 

2 8 0 1 = 1 + 1 

1F(PATBUF < 1 ) ,EQ.409b> PFLAG=I (CHANGE PITCH? 

1 = 1 + 1 

2Vr F 0 TO 1 4 0 

C THIS section ados ?EPUS TO oatbuf phfm called. 

300 [0 31.0 J.LAX=JJ,J8 

r o 305 Jl. A y = 1 ,258 

LAJRUF<JLAX + ( ( 11 A>.-1 ) *25&) ) =2048 

5 0 5 C 0 N T I N 0 R 

51 r continue 

GO TO )bO 
C YES. CUIT 
400 eLDFILL 2 



E NOFILE 3 
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E t\> D r T L L ti 

si or 



f fj 0 
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COEF 




FIGURE A-1 1 
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C. 

C 

C 

r 

c 

r 

c: 

r 

c 

c 

c 

r 

c 

c 

c 

( 

c 



7 I T i FT : COirT.PUR 7^1206 

THIS S ll 3 P 0 U * I N E T^KlS T hi y/ ALuEiS THAT Wt Rf Kf AD BY WINDOW 
A T 0 C A L. C 1 1 L A I F S IMF COEFFICIENTS FOR THE APPROPRIATE 
WINDOW ImaT WAS RfOUrSTf L . IT KflUPMS rurst coefficients 
to r t ,e "Iain program in arfay "a". 

DFFIA n IONS s 

1. A - Ark A v CONTAINING C UEFF I C I LMTS . 

?. f''T y pc - THt valle that corresponds to thl 
WIN pOW OP FILTH dlsirld. 

3 . k - the value ip/T determines the length 

of THE p.-pUl SF I FSPONSL r)F THE WINDOW 

t 

AND thl PUMPER Of COLFKlCltNTS CALCULATED, 

A, L = 2 * K + 2 



SUBROUTINE COEF ( NT Y T E t K , A « L ) 
DIMENSION a ( 1 5 J 

[0 m i = 1,15 

A ( 1 1 = 0.0 

TP (OMTTN'JF 

i ) = 3 . m i 39 

s c a i r = o.o 



L = 2 * Y + 2 

(0 TP ( ] 00 1 2U0 , 300 , t, 00 , 600 . feSO ) , NTYPE 
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0 H A M ! J 1 r, G W 1 M DOW Spi FoTtf 

3f'p ro no i = i,l-i 

Ml) = .5 - # *5*(C0S( (2*Pl*l)/l > ) 

SCA| F = SCALF + All) 

13 p cOMTI^UF 

(-0 TO 50 0 

(. MlNUQW S 1 1 FCTt.C. 

?r P TO pj 0 1 =1 , L- 1 

All) = .54 - ( .4 6* ( ( fS( ( 2*^1* I ) /l ) ) ) 

S 0 A| F = SCALF + All) 

23 P rOMTlMUF 

CO TO bp 0 

<: plackfam windOj sfllc^fd 

6 0 p (0 63 0 1 = 1 ,L"1 

All) = .5*1 CCS( ( ?*P1*1 > /L > ) ) + ( .Ofl*(COSl (4*P1* 

1 3 ) / L ) ) ) 

SCAl F = SCALF + All) 

63 p cOMTIMOF 

pO TO bOC 

F PFCI'aACJI.aR LI Ml) 01*' SFl-f LI F D , 

Hop f 0 410 1 = 1 , L- 1 

All) = l.p 

SCAl F = SCALF + All) 

43r cOMtjfJf 

C CO^P-lilL the scale FAC ION. 

M'O 10 A?p ) -• 1,35 



a( n 



A ( I ) / c 'CAl-t 



bRO COMTI^Dt 

CO TO 7 of' 

c rN akfitraRy set or run f icients way or Typed in. 

600 CALI ASSIGN (5» *KKl : • » -1 ) 

l fl F T A E F I L T c (?i?bo*U«NVAr<) 

AVAR = l 

REAP ( b » N \/ A R ) ( A ( J Y ) , T r = 1»?*K+1) 

F N0F1LE 5 
fcO TO 700 

C PaPKS-MCcLFLLaN COFFFJClff TS / \Rt READ IN FROM FILE, 
fcbp call ASSIGN ( 5 « ’ RFl : 1 1 -1 ) 
hEAn ( b) |\!l- XI , jjPG 
AN r (NplL + l)/2 

RE Ap ( b ) ( A ( I Y ) i 1 Y = j , NN) 

LNDFILE 5 

rO FAC J= 1 , MN 

AfNFIL + 1 -J) = A (J) 

600 cOMTTN'Jf 

7 0 p 



RETURN 
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W IN DOW 




F IGURE A-12 
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C 

( 

c 111! F; : WINDOW. FOP. 741220 

C 

CJ THIS sUBPOUIiNe aILOPS FCl’R POSSlpLt WINDOWS, 

c a parks-mcclf;ll/\n ru ur, or an a'*b. filtf-p to 

c PF SFLLCTFO aMD A VALllF CF K (PROPOKTIONal 10 the: 

F iEnGif OF THp wi N n 0 w) TO F E SE.Lt.CTt P T OR jHt MAIN 

r program to use., 

c 

SU3P0U11NE WTMUOW (I TYPE, K) 
l*' K I T F < 7 ♦ 1 U ) 

10 FORNAl<* • ♦ , TYPF IN NUMBER CORRfSP. TO W I MUOW * . / , 1 X , 

1 .DESIRED: HA m N=1 , I ' f. P- M , = 2 » BLACK, =3, Rt.CT.=n» ,/, 

2 «ARpITKaPy = 5« PARf S- MCCLELLAN = 6 (IT)',/) 

h.EAn (7,?o) f'TYPE 

2 0 F 0 R nr A T ( II ) 

LRITF( 7,30) 

3D F ORMAT (« m'tIIL V A L w T OF K = ? (2K 4- 1= LENGTH) » ) 

FEAn (7,40) ^ 

40 FORMAT (12) 

FLTHRN 
F N D 
C 



C 
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FILIN 




FIGURE A-1 3 
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C 

c Tin r: FILIN. FOR 75U31S 

c 

c 

c this subroutine acts io iririALizt the. filter a t THE' 
C PF.G4 hlMG of THU pROFKAj‘1 OR AT » HF OftSE I p F A VOICED 
L AREA A F 1 F R AF IIMVOILFU AREA. 

F 

C IUITI ALIZA I" ION 

S U 3 R 0 U I 1 1\| F FTL1N ( K , IP i T , S * I S U l w i ) 

i i\i r fger s ( j i+ ) 

CO 30 *c = 1«2+K 
S(<F) r I a b S (1PIT x 2) 

ir continue 

L I SUM = T I < f: LfIxIgTh of- 1 he FILTER. 

IS Jr = M * K * IP11 1RLNFW ISUM, 

RETURN 



FDD 



CALC 




FIGURE A-14 
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c 

c 

r 

r 



TITLE: CALC. TOP 7b(,313 



c This SU-3POU1 WE ic used iith 1HL I* A i n program nonvoi 
C AMD is USED TO c A T C L l A T £ /\ M OUTPUT HQINI. 



C PARAr l Trips : 



c arrays: 



c 

c 

( 

c 

C 

c 

c 

r 

r 

c 



1 . A - CONTAINS THl COEFFICIENTS 

2. DA 1 f<UF . CONI A II S THE INPUT DATA 

2. s - contains the specings peTulem the 

COFFFlCTEiMTS. 

VAR I ARLES : 

J 

1 • I - POlWTfK or 1 HE CURRENT PROCESSING in datbue 
?. k - value that lfTeRmines the no. of coef. 

2. I CUT _ a VAt UE I F TURNED TO THE M/iIN PROGRAM. 

THF. OUTPUT or '|HE FILTH FROM one point. 



SUBROL'IiNe CTC (A,LATUUF.i>,I,K,IOuT) 

n pension a ( 15 ) 

IN I EGER DaIBHF < *T09t > » S < 14 > » OUT 

i r ompl t at ton of in'" output point 

Y = 0.0 IINIT1/U1ZE THL SU V| VAUlE. 

Y = All) * [)ATUUI (1) j P t K T 0 H E FIRST NULl. 

IN DEV = T 

FO ?0 1Y = 4 l ! PERFORM REST OF MULT. 
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iMDrX = I \l U E v - b(]' l -1) 

Y = Y + < A( IV) * DAT i UF ( INOt-X ) ) 

yr (OVITINUE 

lOJT = I F T X (Y + .5) IRUUMfJ THL OUTPUT V^LUt 

M. TURN 

FDD 
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I III L; ADAPI .HOC 

DOf Uvjf'MT/' I TUN | OR NONVOi AMP N0 mVO 2 



ABSTRACT 



rpFSE PROGRAMS JM||E>|ENT IMF ApAPTjVE SYSTEMS 
E UR FPL A T L J 1-M THIS THfSjS. NUNVOl CORRESPONDS To I lit 
ADAPTIVE. OVIKlOao SYS I IF WHILE NOE VO? IS THE ADAPTIVE 
SYSTEM. 

SUpROUTjMES U S p p : 

A. MON VO] 

1. Window (ntypl»k) - this subroutine: allows fop 'he 

USr R TC PRESCRIBE 1 HE VlNllOW FUNCTION, P aHKS_MCCLLLL an 
MITER, OR aN A R RlTRAf<Y Ml TER DESIRED IN TpE PROGRAM A^O TO 
umionatf the VAl.UE OF K f C,P 1 he filter length. 

2. CoFF (l\lTYpr,K, Ail ) - THl SUBROUTINE CALCULATES I HL 



( ;. r r F E" I (_ . T L M T S 


OF THF 


F jl TLR 


PROSCRIBED IN 


the 


WINUOW 


SUPROUTIME. and 


STURPS THI 


M I(J AN 


array "A" that 


Can 


BE OF 



n/>i v Uh nifTNSjosi or is. 



S. CALC ( A , OA'I PUF ,S, I , K, TOUT ) - THIS SUBROUTINE 
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Kt T UKNS 0 MT V A[ Up 0 F MlE OUTPUT PER CALL, jl inPLFriENTS I |-|L 

/'Or PTivr fti tor d 1 ffekt^cl fduaiiun. 



4. FILlM (KiMLWPl,S,ISUl“l) - ThJS SUBROUTINE 

JTTTIaLIZcS THF F t LTF.F Al THE BEGINNING OF UlF PPOGK A M aNQ 
A 1 THE Cl'-SFT Or A v°I cru /'FFA FOLLOWING AN UNyOICtD 
TPF F 1 1 1 E p SpACiM r S APT All SET FQUA|. TO THE NEW pIlCH 
PF RjOl) VAL.OF JMT 1 A I l-Y. 

B. NlnNVU? 



1. WTNDOW (iMTVPEiK) - (SEE E XPL. AN A T l ON A ROVE ) . 

2. COFF (NTYPP ,A, AtL) - (SEE LXPLANaH ON ABOVE). 

3. FTI.IN (KiNrWPIiS«lSlJfl) - (SEE EXPLANATION ABOVE). 

4. r A L c 2 ( 1 . PATBLr * S. 1 1 IOOTi ISTaRT ) - THIS 

SdpROUTlME IS BASIC-ALLY Th£ SAME AS TOf CALc SUBROUTINf WITH 
lor F x (|P 1 T n M that tT PRO V I f ES FOR THI Ovt Rl OAP PROBLEM AND 

coppects Pur it. 

5. Brn ( f start , u . i . f ) - this subkouthe Initialises 

tit 1ST ART ARRa.Y WiiJcH IS USED JN THF L'ETFClION OF 1 HE 



- 217 - 



cM/rpi oat RnopLrM. 



G. RfBLG ( TST/'R I , I , K) - THIS SliPPOljTINE RE- I NI T I aL I L ES 
1 1 r j s t a r 1 array after a plfflr shift has bfEim rfkfokmfu. 

OTi r« PRCSPaMS USED WITH 1|PSL SYST|.-p K S 

(STE THt SHI t Ep , D U C EXPLANATION; ALP PROGRAMS ARE THE S Al^t- . > 

INPUTS TO A P A 1 I V E SYSTEM 

1. STRNAl.DAT - FROM ADUNS 

2. BtTCH. TAB - PRON | J T V 0 2 

output r-Rop i H[ adaptive System 
l. output . dat 

AL PORITPM 

1. THE ALGORITHM IS [ ASICALLY THE SaMe AS THE SHIELuS* 

SYSTEM. DM Y THE EvCt-PlJCfS ARE LISTEP In THIS DESCR IpT IU|\|. 

2. THE PITCH. JAB PIL[ IS INITIALLY STORED IN AN ARKAY 

"P" SO EASIER ACCESS Tf. I H{ PUCp PERIOD VALUES CAN BE 

OPT A INCP . 



6. AH ARRAY "S" TS USED TO HOLD THE SPACING VALUES 

ptTWtEi'i the coefficients. the variable i sum is used to 
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lil NOTE THE I E iM GTH or | Ht f 1 1 F K . 

<+. 1 N MONVD? AN AKK/ T "IStA^T" IS USF D TO HOLD I HE 
J PIT I Ai STAKIjNr WAlJJtS (CR A PARTICULAR filter RASED OH A 

njTfH pfiaon Value, j- t is used jn detecting an overload 

(M’rriTicN. 



b. 


i nr output 


POINT 


IS cO.iPHTtD BY 


the 


difference 


L 0 1 ’ A T I C f> 


USED FOR 


1 HE 


/ fAETIvt FILTLR. 


IN 


N0NV02 1 HE 



uvrPLOAi problem is cuPRfLirn. 

6. THF BlJFFtR IS SHlfTEU IN A SIrIlIAr MANNER To 'HE 
SUTFEDS* PpOGKaN. 

7. DURING UNVOICED StrTIONS THE INPl-'T IS ATTENUATED BY 
IHr VARIABLE, con, TO OHT/jN THE OUTPUT. 

o. when the logical variable, done, is true, ihe 

PPOGRAF HERF01S A SERIES Of CLOSING SlfPS, THEN HALTS. 



t rp 



c 



Till r: N 0 M v f 1 « fur 



750311 



r 

( 

C THIS | R 0 0 R 7 ’“I I'THLF'MLHIS THE NONUN II (IRMLY SPACED 

c adaptive filter without correcting fop ihe ovlklcad 

c PROdLlMS. 

C 

C SUBROUTINES USED: 



I . 


WIND Old/ ( A' T Y Pt ♦ K ) 


c p . 


COLE ( MTYPL « K * A . L ) 


c 3 • 


FILIN ( K,Nc.WP) ,c, ISUM) 


C R . 


CALC ( A,DA1PUF,S« I.K, IOUT) 



c 



c arrays: 




c 1 . 


A - CONTAINS C0| FFIC1ENTS. 


c 2 . 


DAIPUF - INPUT VALUES. 


3 3 . 


OU 1 RUT - TLNHCR/FY OUTPUT ARRAY. 


c R . 


p - used to stcrf pitch periods. 


C h. 


S - USED TO SI Of E SPACJNGS BETWEEN COEL . 



L \! A R 1 A F L. E E : 

L i. DOME - LOGICAL [_.L MOTES END OF SEM FENCE, 

c: ?. ATirM - °EijO I b S UNVOICED AREAS. 



C 3. J - POJNTE-k ir' | ITCH PtRTOD ARRAY "P". 

C 4. TSJf" - LTHGTH Of THE FILTER. 

C 

C II\llTTALI?ATI n N 

UlrNSJON A (15) 

I NTrGE 4 Da I'SHF < 4096) , OUTPU I ( 256) , P< 1000 ) , S( 14 ) 
LOGICAL DONE, ATlfcN 

LATA OAiBUK , OUTPUT ,(-,5/4096 + 0.256*0,300 0*0114*0/ 
L A 1 A 3 VAR , JVAR , KVAR/? i 2, 2/ 
rOYF = .FALSF. 

AT 1FM = .TALFF. 

CALI. ASSIGN < 2 . ’RKljSIGNAL.DAT* ) 

CALt ASSIGN ( 3. ’RKl -.PIT CH. I AP* ) 

CALL ASSIGN (4, *K|<l:OUTPUT.DAT * ) 
rEFlNE FILE •> ( 251 , 256 , U , I V AP ) 

CEFTMf FILE 3 (1000, 2, U, JVAR) 

COM = .3 JATTLNUATjON CONSTANT 
C CHOOSE WINDOW jrpF 

CALI. WINDOW (MTYPEiD 

c calculatf the coefficient? 

CALI COfF (NTYIT.K,/ ,L) 

C. SPECIFY the | F. M G T H OF Tht OOTPUI FIIF. 

( 1 ITT = 251 + l K*2 ) + 1 

I EFTNfc. FILE <> (LIMIT ,256, U«KVAP) 

C FILL THE PITCH ARRAY. 



ro 300 1A 



i ,1000 
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\ L Ap ( 3» J\/AR,EIMP = H 0) P ( J A ) 


iro 


CONTINUE 


1 3 r 


( OMTTN'Jl 

o = l. ! J IS THL HOllIEK IN THE PITCH PLROp AHKAy, 


C THIS 


PART DETER K i pFS IT L F HAVE A SILENT AREA AT 


L THt. ( 


F GINN I Ng Op THL F 1 1 L . 


1?P 


IF (P(J).GF.P) CO Tc 200 ! SILENT AREA? 

F = I A US ( P ( J ) * 2) IVES, SET THF OUTPUT = 0. 
CO 150 KO = 1 , ( N-l ) ,? 

KA = KD 


1?0 


IF (KA.LE.25S) CO 1 c 140 
KA = KA - p5s 
(-0 TO 130 


1 M p 


CUT PUT (KA) = 2048 
OUTPUT ( K A + 1 ) =0 

IF (kA.nE. 25S)G0 TO 150 !IS IT T I ME TO WRI1F A pLC.? 
LR1TF ( «H KVAR) < OUT PUT ( KB ) , KB=1 , ?56 ) !YFS. 


1 bp 


continue 


C PUMP 


the counter. 

o = U + 1 




GO TO 1?0 



C PITCI PEP 100 W a S NOT A b ] l ENT 0*'E. 
2f’0 IF(P(J) ,GF»4P0) GO 1 0 3000 
( VOICE l, Sri yp fUF FTLTfl SPACING. 
ASSIGN 350 TO KXX 



CAL! FILIN (K,P( J) , S , I SUM) 
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i rt asure'tfnts section 

IR r: lStlM /2^0 

}F ((IR * *56 ) .ML. 1 Sl |p/ 1 ) IR = IR + 1 
( SET (j( J FOR 1MITJAI i\ F A F 1 1\ f INTO BUFFFR. 

IVAP = KVAR - IR 
CO TOO KE r 1,16 

F t A P ( ? * I V A R ) ( u A I L>Uf ( KF + ( ( Kfc -1 ) +256 ) ) t KF r 1,?56) 

p f 0 \IT T M J£ 

( SET up IMF 3pFFE R POINTER 

IF (KA.EQ.^5^) KA = -1 ! RESF1 THE POINTER 

! TO THE LAST OUTPUT , 

I = (IK * 256 ) + ( Ka +• 2 ) 

c ro T | iF OtiTPJT CALruLATIOfS 
55p CALL CAlC ( A , OATBUF « S * 1 « K « I OUT ) 

C COMPUTE THE ADDRESS FUR THE OUTPUT POINI. 

56 P HMD = 1 

5 6 f. IF (IN0.LE.2S5) 00 TO 570 ! DO WF HAVE A GOOD AUDREYS 

I ND = IpD - 250 ! MO Go RECHECK. 

(-0 TO 3 65 

57 n 1 F ( A T 1 E K ) GO TO 575 ! I S THIS A VOICED SEGMENT? 

C U 1 P U T ( I N D ) •= I OUT | YES, GET CAl CULATEU VALUE. 

FC TO 5p 0 

C ATTLMATF the INPUT. 

5 75 GUI PUT ( I (JO ) = ( COM* ( (■■/ TBUF(I)-?046))+2048 

56 0 • CUTpi.iT ( lNO+1 ) = 0 ! SECOND CHANNEL r 0. 

f IS IT TICE IQ r'Uvjp f HE DU PUT BJFFt R? 
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IP ( INO.NE* 2 C B) GO 10 400 

UR1TF( 4*KVAR,F(4D = ?0g 0 ) ( OU I PUT ( KP ) , KR = 1 , 256 ) !YpS. 

400 1 = 1+1 

0 IS Tin RE A PITCH r/\ KK? 

IF ( D A I F U F ( I ) . L T . 4 0 7 F ) GO 10 12?0 
bf'O REWpi = R(J) * 2 ! Y t S , GET A Mb UJ VALUE OF PITCH. 

^ = J + 1 ! BUMP Hi FITCH PERIOF POINTER. 

C /'RE V F EMTLK1NR A SILENT AREA? 

] F (NEWPI.t T.O) GO 10 1600 

c ro, are we entering ah ui. voiced area? 

IF ( MEWPI.GT.4U0) be TO 1BO0 

c: mo« iiERrroKf we rRL trip ring a voiced arfa. 
c chance the filter. 

ASSIGN 35n T 0 KXX 

C IS Tlif FTtlEp BEING INITIALIZED AFT l R AN UNVOICED AREA? 

if (atteN) gr to eoo 

CO TO 900 ! THE FU T[R IS NUT BEING RE- 1 N I T 1 flL I Z ED . 
BOO AllFP = .FALSE. ! Cl FAR THt FLAG. 

CALL FIllM ( v . NEWP 1 , S * 1 SUM ) 

(0 TO 1?20 
9”G R A A = 2 * y 

3 000 ][ (KAA.FD.i ) GO 10 1100 

b ( AAA) = SCKAA - 3) ! SHIFT TFIE SPACIfGS 

F A A = A A A -1 

• go to ’ron 

lirn s ( 1 ) = NEWPI ! SET Iff NEW PITCH PERIOD SPACING 
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]SJr = o IPO^rUTL T|F L E f\j G 1 H OF THE FILTER 
LO 1200 y i = 1 , ? * K 

is Jr. = jsum + s(Ki) 


3 20 0 


C OVITINUE 


1??0 


I = T + 1 



( JS IT TlMF 10 SHIFT Tut [ OFFER? 





IF ( I . LT . 4U9M GO 1c 140 0 


0 Yf S, 


shift fhe puffep. 
ir = i sum / ?56 

IF (( 1**256) .Nt.ISUr) IR = JR + 1 
IVAR = I V A R - IR 
10 1300 K? = 1,16 
IF ( IVAR.GE.P52) ro TO 2900 

FiLAn (P'lVAR) (PATFHFFK3 +- ( ( K2-1 ) *25f ) ) * K3 = 1,25b) 


3 3(1 p 


CONTINUE 



t RFSEl 1 



3 3.1 0 


1 = (IK * ?5F, ) + 3 


l‘i 0 0 


CO TO KXX , (350,360,1840) 


0 THIS 


AREA RESFRVcr FOR SilENT AREA LATER. 


IF 0 0 


RE /PI = IABS(NEVRJ) ! GET RID Of MINUS S I GN 
covf = .true. i se t done flag. 

ro to icon 



0 UM/OIftJ ARE A 



l v r c 


1 = 1+1 sblfp the pouter. 

IF (F'E/pI.l. E. ISDN) ( 0 10 2000 


3. 8 ? 0 


F 7 = 1 ! St T UP A COUNTER. 
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1B?5 K A 6 r 0 

i: are wf through shtfting? 

1K30 JF (KAd.LO * ( ( 2 * K > - h 7 ) )GO TO 1835 
KNEW = 2 * K ~ KAP 
l< 0 -P = K N r w - 1 

S(KMrW) = S(KOLP) ! no, CONTINUE 

HAB = KAB + 1 
ro to iP3n 
1835 Kti = 1 

I837 ASSJGN 1040 to KXX - 

CO TO 350 ! GO ComUlf AN OUTPUT, 

18 40 Kfi n KB + 2 

IP (K 6 .LT.S(J ) ) GC IP 1837 ! AKE WE THRU? 

K 7 = K7 + 1 | YES , BUMP K7 

c enter the attfnuation procedure 

1B5Q A T 1 r r .! = , TRUr. !SET THE ATTpNUATlUN FLAG. 

IF { pONp ) GO TO ISOC ! AKE Wt ALMOST THROUGH? 

ASSIGN 360 TP KXX ! |. r , THIS TS ONLY fiN UNVOICED AREA. 
CO TO 5f0 
3 V 0 0 3 NO = I 

0 ro W( HAVE A SUM.^PLE IN|_? 

1V?0 IF (IN0.LE.255) GO IC 1950 

] NO = 1ND - “55b ! NO « TRY AGAIN. 

CO TO 1920 

1980 . CUTrUT(lNp) = <i048 (SILENT ARtA, THEREFORE 



! SL 1 OUTPU1 



2048 (OR 0) 
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fUTptlT (IMP i 1 ) = (. 

If ( IND.NE«255) DO 10 1970 ITIPT 10 DUOP OUTPUT BUG? 
URlTrl “M K\/AR,flMn = 2,00 0 ) < OUTPUT ( KB >♦ KP=1 , 2?6 ) 

ifJO = -1 ! K[ SE 1 POINTER 

1970 I NO = IMP + ? 

GO TO l92n ! C( FT . UNI II. IMP OF FILE. 

C THIS Ah£.A Will. BE PlALIiMG WITH UNVOICED SEGMENTS WITH 
C TSUM>M. WPI . 

2 0 0 0 GO TO 3000 

2900 TO ?91 0 I LAX = K2,1l 
TO 2906 JLAX = 1*25( 

LAIRIJFIJLAX + ( ( 1LAX-1 ) *25b) ) = 2048 
29n5 CONTINUE 
2910 CONTINUE 

GO TO 1310 
3000 E NOFILE 2 

f NOFILE 3 
FNDFILE 4 
STOP 

END 



NON V02 




FIGURE A-15 
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N0NV02 




© © 



FIGURE A -15 (CON'T ) 
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N0NV02 




F IGURE A -15 (CON'T.) 
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c 

c 

c 

r 

r 

L 

C 

.C 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c: 

c 

c 

c 

c 

c 

c 

c 



TITic: NoiMVOp.rUK 75U3PA 



TI'IS \ R Or R A ^ IVIHLrMLN'IS T HE NONUINJIRORMLY SPACED 
ADAPTIVE FILTER WTTH the CORRECUON for the overload 
PROIJlf tf . 

SUPROLTIMES used: 

1. WJMpOW ( rT YPE * K ) 

P. COEF ( MTYPLt Kt A, I ) 

3. FILIN ( k, NLPPI t St ISUM) 

4. PEG (PSTARTtSfliM 

5. CALC 2 ( At DATHUF f St 1 1 K« 101'Tt 1ST ART) 

G. REBEG ( Id/vR I t 1,K) 

aprays : 

1. A - CONTAINS COEFFICIENTS. 

p. p - used for rescaled loff. during overload. 
3. DAVpUF - CONTAINS 1NPUI. 

H . OUTPUT - TlEPOK/PY OUTPUT ARRAY. 

E • P - CONTAINS f 1 i t r H PERIODS. 

f. S • CONaTInS SP/ C -TNG FOR FILTER. 



C V AR i AF E E S : 
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C 



]. no\ir - logical rrNOTEs f-: m n of sfnifncl 



2. ft T T [ N - LOGICAL | EN0TES IJNVf) I LED AREAS 



c 



16. J - POliMTF'K IN HE PITCH I’tRIOL) ARRAY "P" 



C 



A . I S Ll F _ LENGTH 0^ THE FILTER 



C 



C 



initialization 



L TENSION A(15 )i B< 15) 

JNITOEK DA TBI'F( 4 096) , OUTpUl ( ?56 ) , P ( 1 0 0 0 ) , S ( 1 4 ) 
INTEGER IRIAOT (15) 

LOGICAL done, aTIEN 

FATA DATBUE, 0U1PUT,L ,S/409b*0, 256*0, 1000*0,14*0 / 
data ISjArt, 0/15*0,15*0.0/ 

L A 1 A I VAR, JL/AR i KVAK/P, 2, 2/ 

CATa OONF, ATTLN/.FaI SE. , .FAL SE./ 

CALL ASSIGN ( 2 , ' RK1 ; SIGNAL. OAT • ) 

CALI ASSIGN ( 3 , ’KKljPlTCH. IAB’ ) 

(.ALL ASSIGN (4» 'HK1 ICUTPUT. DAT’ ) 

LEFTf'E FILE 2 ( 251 , £ 56 , U , I V AR ) 

I.FFINF FILE 3 (1 U00 , 2 t U, JVAR ) 

COM = .3 

C riiOOSF /JTNDOB t t pr 

CALL W1NDON (MTYPLth) 
c calculate the cuEFFicirri s 

CALI COfiF ( N T Y I J E » K » ^ , L ) 

C SPECIFY THE LENGTH OF THl outplji Flip. 



L 1 Y T T 



2 5* + ( K* 2 ) + 1 
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[ LFTK'fc FILL <! (LJril ,2b6.UtKVAK) 
r fill the. pi i ch array . 

LO 100 1 A = 1 ,1000 
HL An ( 3 * JVA R , E'HjO = HO) P(1A) 
inn l 0 \! T T N U f 

110 cOMTTNUr 

t 

o = l i J IS 1 HL POINTER IN THE PITCH PER OR ARRAY. 

C THIS FART OETERfipTS IF t.F HAVE A SllFNT ARF A A1 
C THE [fGlNMI.MG Of THE ULl, 

IPO IF (P(J).GE.f') 00 Tc 20.0 ISILEMT AREA? 

f = I A c3 S ( P ( J ) * 2) ! YES, SfT THE OUTPUT = q. 

L'O 150 hD = * , ( N-l ) , ? 

KA = K 0 

130 IF (KA.LF.255) OU To 140 

KA = K A - ?5F 
00 TO 130 

1^0 CUTOUT (KA) r 2040 

OUTPUT ( KA + 1 ) = 0 

IF (kA.kE.255) 60 To 150 ! I S II TINE TO WRITE A RFC.? 

UUTF< 4 • KV/'R) (OUTFLT(KB) .*B=1 .?56>) |YES. 

150 COMTINJf 

( PUNP THE COJMEK. 
w = J + 1 
(0 TO 1 2 0 

r PITCH PERTOO Was uor A S i l ENT ONE. 

20 0 IF ( n ( J ) . C E • 4 no ) 00 -) C 30 0 0 
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< VOICfL t STl UP I nr FILTL| SPACING. 

ASSIGN 350 TP KXX 

CALL F1LIM (K,P( J) 'S' I SUM) 

L PEASUhC’irPT S SECTION 
JR = I SUM /2 r '6 

JL ((IK * 256) .ML. 1*1 M) IR = JR + 1 
C St T UL FOR INITIAL REAP I N C INTO BUFFER • 

I V Ap = KVAK - IP 

to 300 Kt = 1 , IF 

RLAP(?’IVAK) (UATPU|. (KF. + ( (l<r-J >*256) ) ,KF = 1*256) 
30 p cOMtINUf 

C SET Ur IMF But Ft. R PU INTEL. 

IF (KA.FG.25s) KA = -1 j RLSFl THE P 0 T N T F R 

! OF THE LAST OUTPUT. 

J = (IK * 25M + ( K„ + 2) 

C TMTI/-L17P I HE /START ARRAY 
CALI BEG ( LSTAKT.S' j ,K) 

r re thl output calpulat iors 

3 b 0 CALL CALC? ( A , 0 A1 PUf , S . I , K * I OUT , I ST ART ) 

C COMPUTE THE ADDRESS FUR HE OUTPUT POIN1. 

36 r IND = I 

36s JR ( IND.LF.2 r 'f>) 60 T p 370 ! DP U'F HA VF A GOOD ADDRESS 
1 IM □ = J ND - ?56 !N0 60 RECHECK. 

CO TO 3p5 

3/p iF(ATlEt) GO TO 37b ! T S THIS A VOICED SEGMENT? 



CUIPHTLIIJD) - I OUT 



!YES, GFT CALCULATED VAlUE 
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GO TC 3f,0 

C ATTLI l A 1 r 1 HE. INPUT. 

37f CL) 1 PUT ( I r J 0 ) = c cor =4 (DATBUFC I)-2046) ) + 2048 

38 p GUI PUT < jhD + 1 ) = 0 I SECOND CHANNEL = 0. 

C TS IT Tlr-’F TO nU V|P I ML OUTPUT QJFFEK? 

IF (IM0.ME.P8b) 00 10 400 

LK1TF( 4 ’ KVAk, El'jn = 3l'u0 ) ( OU I PUT ( KB ) * KB = 1 . 256 ) ! YES. 

4 II P 1 = 1 + 1 

c is there a hitch nark? 

IF ( P A T B U F ( I ) , L T » 4 0 7 b ) GO 10 1220 
blip kEWPl = P(j) * 2 ! Y [ S , GET A NEW VALUE OF PITCH. 

v- = J + 3 JBIIMH THE rilCH PERIOD i j ointlr, 

C ARE he entering a SILENT area? 

)F (NEWFI.LT.O) GO i 0 1600 
c no. are we entering am ur voiced area? 

IF (MEWPI.GT.4U0) be TO 1900 
r NO, IFEREFOKt WE "Re t NTf PING A VOICED AREA. 

C CHANGE the FILTER. 

assign 35n TO KXX 

c ts t i f fjlter being initialized after an unvoiced arfa? 
IF ( A T 1 f N ) GO ro 80 G 

GO TO ypo !1HF FllltP IS NOT BEING Rf- I NI T I AL I ZfD . 
800 ATlTN = . F«L^E . E C i t AK THt FLAG. 

C TINE TO SLIDE THE 0UFFFK? 

IF ( I . Gl . ( iMEwpI*i? + K) ) GO TO 820 

,J )0| D = I ! Y i S , SAVE THE POINTER. 



- 235 - 



1 H A P = HEWP1 * 2 * p 
1 RA = ( ] It A K - 1 ) / 2 bt 

IF ( ( ,KF. ( Jf AP-I) ) IRA = IRA + 1 

IV A P = ] V A R - LG - 1 P A 

CO 810 JOY = 1 1 10 

RE.An ( ? • I V A R ) ( DATl'LF ( 1DX + ( ( lflY-l ) *2bfe ) ) , TDX = 1 « 2S6 ) 

Ain (OUTINUf 

C' RTSET Trip POINTER 

] = JH0| D + <200 * J R A ) 

020 (ALL FILIM ( Ki (NFWH/2), S,TSUR!) 

CALL BEG ( ISTAKTi S» j , K) 

CO TO 1220 
VOp K A A r 2 * K 

.1 000 IF(KAA.FO.I) GO 10 j.100 

? ( K A A ) = S ( K A A - 1) ISHIFT THL SRACINGS 

K A A = k a a -J 
GO TO 1000 

nno ? ( i ) = newpi fsri iir new pitch period spacing 

iSJP! = 0 ICOI^PUTL T j F LEMG I H OF THE FILTER 
CO 1200 Kl = 1 » 2* ^ 

J S Jr = jSlJ'W .+ S ( Kl ) 

1200 C 0 M T T M 0 [ 

CALI BEG ( I START « S * } + 1 « K ) 

1220 1=1+1 

C IS IT TIFT 10 SHITT THf [OFFER? 



IF < I » L J ,4098) GO 1c 1400 
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c yets * shift ihe pjrFtP. 

IR r I SUM / P56 

IF ( ( 1R*?56) . ML. 1SU| ) IR = IR 4 1 
]VAr = i var - ir 

r 0 3 3 0 0 K 2 = 1,16 
IF (lVAR.Gt.P52) GO TO 2900 

I < E An (2*IVAR) ( 0 A T b (_ F ( K 3 4 ( (K?-l )*25&> ) «K3 = 1«256) 
1300 (OVTIMUf 
G prstl I 

1310 1 = ( I K * ?56 ) + ] 

C PFSEl T Hr ISTARI ARRAY. 

CALL RtpEG (TSIAHT,},K) 

1M00 GO TO KXX, (35Ci,3GO,3BMO) 

0 THIS AREA RESFRv/EP FOR SiLENT AREA EaTCR. 

1600 rEa/pI = I APS ( E.tWPl ) ! Gt T Rip OF M JNUS SIGN 

r c Mr = .true. istT done flag. 

(-0 TO lpoo 
G UP VO 1 c E 0 area 

1=1+1 ILbFP THE POINTER. 

IF (NlWpl.l E. I SUN) LG TO 20 00 
1 6? 0 h 7 = 1 !S[T UP ft COUNTER. 

1(>?5 h A3 = 0 

g ape u through shifting? 

1 A 7 0 IF (KA3.E0. ( ( ? * K ) - K 7 ) ) GO 10 1(3?5 

h N E 1-1 = 2 * K - KAP 



hCLO = KNEW 



- 237 - 



M<f!FW) = S(KCLP) ! NO, CONTINUE 

h A 3 = K/'R + 1 
fO TO lG3n 
1R3E iSJM = fi 

ro 1.836 1 Tf r i,2*K+1 
] S JN = ] SUP" + S < I T L ) 

10?£. fOMTTNUf 

TALI BtO ( 1STAKT»S» i+l« K) 

Ke = i 

]<■? 7 ASS TOM iB-4 0 TO KXX 

00 TO 3f,0 ! GO COMPUlF AIM OUTPUI. 

I 8 U 0 K 6 = K 3 + 2 

IF (K8.LT.SM)) GO i 0 1837 ! ARE WE THRU? 

K 7 = K 7 + i ! Y l S i BUMR K7 

IF ( K7 - ( 2 * K ) ) 1825,1835,1850 ! ARE wE THRU? 

c fmter tht aitenuatjom prccluure 

1850 ATTFM = .TRUE. ! Sf T 1 HE ATTENUATION FLAG. 

IF (DOME) GO TO 190c ! ARE Wt. ALMOST THROUGH? 

ASSTCIM 30n TO KXX !((, THIS IS ONLY aN UNVOICED AREA. 
GO TO 300 

19po 1 N 0 = 1 

c pr wf have a suitable n«i ? 

1920 IF ( I N 0 , l. E » 2 5 5 ) GO ',0 1950 

] NO = J r D - ?5t> !M0, TRY AGAjM. 

GO TO 1 92 0 



195 0 



GUI PUl ( j N 0 ) 



2 04 8 



fSILFMI AREA, I HERLFORL 
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! St. I OUTPUT = 2048 (OR 0 ) . 

rUI PUT ( I MO + 1 ) = L 

JP ( I NO, N'T. 2*5) 00 If 1970 ! T 1 f F 10 pUMP OUTPUT BUP? 

U«ITF( 4 tK\/AR,FND = 3000 ) ( OUT PUT ( KB ) , KB = 1 , 256 ) 

TUO = -] ! K L S E 1 P01MTFK 

1070 1ND = 1 |\;C + ? 

TO TO 1 52 0 ! C C fv T . UNIIL tND OF FILF. 

C THIS aH E/*. FILL be DlA.LINc WITH UNVOifEU SEF^EHTS WITH 
F ISUFX' L A. PT . 

2000 W [< I T E ( 7,2200 ) 

2200 FORMAT (' » ,‘PROUFAF EXITEU FROM 2000*) 

C-C TO 3000 

29C0 TO 2910 I|_Ax = K2,lfc. 

LO 2905 JLAx = l«25t 

C A 1 RIJF ( JL A * + ( ( 1 LAX -1) *255) ) = 2048 

2905 CONTINUE 
2910 CONTINUE 

C-C TO 1310 
5000 F N D r I L. E 2 
FNDFTLE 3 
[ NDFILE 4 
9 TOP 
l f D 



BEG 




FIGURE A-1 6 
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C 

C lUir: BEG.fOK 7bo323 

C 

C 

C THIS pR0gRA |W I IS a SUpKOUlTWL FO* T H[ NONyOZ.FOK PROGRAM, 
C tt 1 1 T 1 1 A L I 4" F S I HP I S I A K T ARRAY WHEN CALLF'O. 

c 

C INITIALIZATION 

SU3poU1 i ME BOG ( PST/ f'T i S, I » K ) 

IN«EOF K PS I APT ( 1 b ) . S( 14 ) 
c 7 FRO THE PSIaRT /\kr,\Y 

ro 1(1 iu = 1 , 2 + K + 1 
FS! APT ( IU) = 0 
IP fOMTTNJF 

C I N I T I /, L 1 7 F r H E ARPAT WI1|. THE SKATINGS. 

KND = I 

PS 1 ART (1) = I 

LC ?0 J r ^ , ( 2+K ) + I 

hND r KND - S ( J - 1) 
f-START (J) - KNO 

70 CONTTNOf 

hEIilRN 
fc NO 
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REBEG 




F IGURE A-1 7 
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( 

f 

C llll.r: KEPFP.HPK 7503P6 

c 

c TITS SU3P0U r I Nr RF-1H1T1/ L IZES I HE JSTaKT AHRAY AFTER 
C A shirt If 1 PE BUFFlR H/<. BEER) PERFORMED, 
c 

c initialisation 

S H 3 R 0 U 1 INF REPLF ( IsTAKT , J » K ) 

IMFF-ER ISI'ARTUb) 

C PFRFORF THE [ X CHa^'GL OF UE ELEMENTS IN ISTART. 

ro ?o J = i, ?*« + i 

ISTART (J) = I - U,096 - lSTAKT(J)) 

?o comtinof 

RETURN 

fnd 
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CALC2 




FIGURE A-1 8 



- 244 - 



C 

(' 

( 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

(• 

r 



tiiif:: calc?. for 



7 5 R 3 ? 4 



THIS SUBROUTINE IS USED h ITU THL PAIN PROGRAM NONVO? 

afo is used to caiculatl an output poini, it also 
rfTLRMlMrS If' THE COEF E 1 <. I EN TS NFFfT TO BE "TURNED OFF", 
r A R A F ft lrs : 

arrays: 

1. A - CONTAINS thl coefficients 
?. OATbUF - CONTAINS 1 HE INPUT DATA 

3. S - CONTAINS Tht SPACIN 6 S BETWEEN * HE 

COF FFJCTENTS. 

4. TSIaRT - CONTAINS THE INITIAL LOCATIONS 

OF THE rOLF. f CP THE PRtSFMT FILTER. 

VARIABLES : 

i. I - pointer of ife current processing in datbue 
?. k - value that letermnes the no. of coef. 

3 . TOUT - A VALUE KETURNEU TO THE MAIN PROGRAM, 

THE OUTPUT UF TEC FlLIER L ROM ONE POINT. 



SUBROUTINE CALC? < f , F A TBUF i S . J < K * 1 OUT . 1 ST ART ) 
l I VirNSloN A ( 1 5) i B ( 3 *=• ) 

INTEGER CMSUF (40^>f ) »S(14) , TSTARK15) 

LOGICAL SC ALT 
• S'CAl E = .FALRE. 

L SEE IK A A* Y BAD SAMPLES h/.VE BEEN REACHED. 
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^l\iD = J 
H ( 1 ) = A ( 1 ) 

LO 10 Ip = ?, (2 *Y) 4 1 
JMD = JND - S( iP -1 ) 

C rc WF H A v r A COt-FFILTtUT THAT HAS HASS L U JTS 
t: AREAS? 

II ( JMU.GE. ISThPT( If -1) ) GO TO 0 
F (IP) = A ( I P ) 

GO TO 1 C 

o hup) = o.o ! v t. s i turn off the coef. 

sen r = . TRur . 

30 CONTINUE 

IF (.MOT.SrAlF) GO 10 18 !UO WE NEED TO RESCALE? 
SCA| F = .FALSE. IlfS, RESCALE, 
sc al = r.n 

0 COMPUTE THE RESCAI F FACl( P 
ro 15 IT = 1 , (2*K) + 1 
PC Ac = SCAl + b(H) 

15 (.ONTTN'Jf 

C RESCALE THE COEt F T E C T EF T s 

(0 17 IU = 1 , ( 2 *K) 4 1 

F ( IH) = B( I U) / S C A L 

17 CONTINUE 

0 CCMPUTAlTnN OF T|ir OUTPUT point 

18 y = 0 . 0 

Y = P ( 1 ) * DAT bUF ( 1 ) 



indfx = I 

LO ?0 ly = ?, (2*K) + 1 
INDFX = IM|.'EX ~ *>(IY -1) 

Y = Y + ( P ( TV)* UA1 [ LF( INDtX) ) 
CONTINUE 

I OUT = IFTX (Y + .5) 



RETURN 
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